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PREFACE 


This  final  report  covers  the  work  performed  during  fiscal  year 
1975  in  applying  and  further  developing  the  Transportation  Systems 
Center  (TSC)  Instrument  Landing  System  (ILS)  performance  predic¬ 
tion  model.  In  addition  to  the  model  applications  and  development 
reported  here,  the  TSC  group  also  undertook  a  limited  scale  model¬ 
ing  effort.  The  results  of  this  scale  modeling  effort  will  be 
released  in  two  special  reports  to  be  issued  early  in  fiscal  year 
1974. 

The  work  reported  on  here  was  performed  by  members  of  the 
Modeling  and  Analysis  Section  of  the  Electromagnetic  Technology 
Division,  Transportation  Systems  Center  for  the  Category  I  and  II 
ILS  Section  of  the  Terminal  Navigation  Branch  of  the  Systems  Re¬ 
search  and  Development  Service,  Federal  Aviation  Administration 
CFAA). 
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1.  INTRODUCTION 


The  objective  of  the  Instrument  Landing  System  (ILS)  Perform¬ 
ance  Prediction  Program  is  to  provide  a  relatively  inexpensive 
scientific  means  by  which  ILS  performance  may  be  predicted-  The 
prediction  is  accomplished  through  the  use  of  a  physics  model 
based  on  electromagnetic  scattering  theory  whose  equations  have 
been  programmed  for  use  on  a  high  speed  computer.  The  model 
predicts  such  things  as  the  ILS  performance  that  would  result  from 
a  proposed  addition  at  an  airport  of  hangars,  terminal  buildings 
or  control  towers.  Suggested  modifications  to  the  proposed  build¬ 
ing  plan  may  be  made  following  the  prediction  study.  The  model  is 
also  used  to  predict  comparative  ILS  antenna  array  performance. 

This  prediction  is  useful  if  one  is  in  doubt  as  to  the  ILS  needed 
for  a  new  runway  instrumentation,  or  for  tne  upgrading  to  more 
stringent  FAA  requirements  of  an  already  instrumented  runway,  or 
is  useful  as  a  means  of  sti  ’ving  a  problem  site  which  has  not 
produced  an  acceptable  course  structure  with  the  existing  ILS 
antenna  array.  The  need  for  expensive  testing  of  different  antenna 
systems  for  such  sites  is  then  considerably  reduced  with  the  use 
of  the  ILS  performance  prediction  model. 

The  users  of  the  ILS  Performance  Prediction  Model  will  be 
personnel  from  the  Federal  Aviation  Administration  (FAA).  The 
model  Jill  serve  as  an  aid  to  the  FAA  personnel  in  helping  with 
the  prediction  of  ILS  performance  when  modification?  to  the  exist¬ 
ing  airport  environment  have  been  proposed.  The  FAA  should  then 
be  in  a  better  position  to  recommend  any  necessary  changes  as 
indicated  by  the  model. 

A  Phase  I  ILS  localizer  computer  program  package  for  localizer 
performance  prediction  has  been  developed  and  has  been  delivered 
to  the  FAA-  This  program  is  capable  of  predicting  localizer 
performance  in  a  medium  complex  airport  environment  where  the 
scatterers  can  be  represented  as  large  rectangular  or  cyiindircsi 
perfect  reflectors.  These  generally  are  not  important  limitations 
as  most  important  deregators  are  typically  the  large  nearly 
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perfectly  conducting  metallic  or  metal  reinforced  concrete  hangars 
and  terminal  buildings.  For  structures  which  are  not  perfectly 
reflecting,  the  model  gives  a  worst  case  result.  The  basic 
electromagnetic  scattering  theory  of  this  Phase  1  local iter  per¬ 
formance  prediction  model  is  presented  in  this  reports 

The  computer  program  written  for  this  physics  model  has  been 
structured  to  allow  easy  variation  of  key  parameters.  For  example, 
the  performance  of  a  particular  localizer  antenna  array  in  a 
changing  airport  environment  say  be  predicted  in  order  to  determine 
the  effect  on  the  localizer  performance  of  different  building 
locations  and  orientations;  and  the  relative  performance  of  dif¬ 
ferent  ILS  localizer  antenna  arrays  may  be  predicted  for  a  given 
airport  environment.  This  parametrization  of  the  model  allows  a 
variation  of  either  the  antenna  system  or  the  airport  environment 
for  systematic  comparative  study  predictions. 

A  number  of  such  studies  were  made  this  year  with  the  Phase  1 
ILS  localizer  performance  prediction  computer  program.  These 
applications  of  the  model  to  specific  airports  are  discussed  in 
Section  5  of  this  report. 

This  report  is  essentially  divided  into  two  parts.  These  are 
an  applications  part,  mentioned  above,  and  a  theoretical  D3rt  in 
which  localizer,  glide  slope  terrain  and  parabolic  antenna  scat¬ 
tering  theory  are  developed. 

In  Section  2  the  physics  basic  to  the  development  of  the  ILS 
localizer  performance  prediction  model  is  developed  with  refer¬ 
ence  to  an  earlier  Transportation  Systems  Center  fTSC)  publication. 3 
This  earlier  work  is  further  developed  to  allow  a  determination  of 
localizer  performance  when  the  localizer  radiation  is  scattered 
from  a  rectangular  perfectly-conducting  slab  with  arbitrary  tilt 
angles  and  heights  above  the  ground  plane.  In  the  same  section, 
the  localizer  scattering  theory  for  vertical  cylinders  (simulating 
water  towers  and  other  cylindrical  scattercrsj  is  developed. 
Appendixes  A.  B  and  £  further  extend  these  results-  Appenoix  A 
develops  the  theory  of  scattering  Fro s?  cylirlrical  objects  raised 
above  the  ground  plane.  In  Appendix  B  the  scattering  theory  for 


triangular  shaped  objects  with  arbitrary  tilt  angles  and  heights 
above  the  ground  plane  is  presented.  This  latter  is  useful  in 
shortening  the  running  time  of  the  computer  program  particularly 
when  complex  shaped  buildings  which  can  be  broken  up  into  rectangu¬ 
lar  and  triangular  shaped  parts  are  present.  A  method  of  predicting 
localizer  performance  when  there  is  a  hump  in  the  runway  is  devel¬ 
oped  in  Appendix  C.  Since  this  is  a  relatively  coaaon  occurrence, 
its  addition  to  the  Performance  Prediction  Model  extends  the  appli¬ 
cability  of  the  aodel. 

Section  5  presents  the  theory  of  1LS  glide  slope  terrain 
scattering  under  conditions  of  important  large  scale  terrain 
variations.  Applications  illustrating  this  glide  slope  terrain 
scattering  theory  are  presented. 

Section  4  presents  the  method  and  theory  used  to  analyze  the 
Texas  Instruments  parabolic  localizer  antenna  for  inclusion  in 
the  ILS  localizer  performance  prediction  models  Appendix  E 
contains  a  short  discussion  on  localizer  array  patterns  when 
mutual  coupling  is  present. 

Section  5  discusses  the  applications  of  the  ILS  localizer 
performance  prediction  model  to  different  airports  studied  this 
year.  Finally,  in  Appendix  B,  glide  slope  fuselage  scattering  is 
discussed  and  a  method  of  evaluating  the  scattering  functions  is 
presented. 


2.  LOCALIZER  SCATTERING 


2.1  INTRODUCTION 

The  techniques  for  investigating  ILS  multipath  problems  which 
were  developed  in  the  TSC-FAA  report  "Instfunent  Landing  System 
Scattering*’1  by  6.  Chin,  L.  Jordan,  D.  Kahn  and  S-  Morin  have  been 
applied  to  two  previously  untre-'ted  scatterers.  Specifically,  we 
have  developed  and  incorporated  into  our  basic  localizer  coaputer 
program  Xoraulas  describing  the  scattering  characteristics  of  an 
elevated,  tilted  rectangular  slab  and  a  right,  circular  cylinder. 
The  detailed  calculations  of  these  scattering  formulas  are  pre¬ 
sented  in  the  following  two  sections. 

2.2  LOCALIZER  SIGNAL  SCATTERING  BY  .%S  ELEVATED  RECTANGULAR  SLAB 

In  Reference  1,  localizer  sigi al  scattering,  by  a  flat, 
vertical,  rectangular  wall  resting  on  the  ground  was  investigated. 
Ke  new  examine  the  more  general  problem  of  localizer  scattering 
by  a  flat,  rectangular  slab  which  is  elevated  above  the  ground  and 
tilted  away  from  the  vertical.  The  base  of  the  slab  is  assumed  to 
be  parallel  to  the  ground  plane.  It  is  further  assumed  that  the 
slab  is  perfectly  conducting. 

The  geometry  of  the  problem  is  illustrated  in  Figure  2-1. 

The  x-y  plane  is  the  ground  plane.  For  simplicity,  it  is  assumed 
that  the  ground  is  perfectly  flat  and  perfectly  conducting.  The 
x-axis  represents  the  center  line  of  the  runway  while  the  z-axis 
Is  perpendicular  to  the  ground-  The  unit  vectors  in  the  x,  y,  and 
z  directions  are  denoted  by  ey,  and  e_  respectively.  The  unit 
vector  n  is  parallel  to  the  base  of  the  slab  while  €  is  perpendi¬ 
cular  to  n  and  lies  in  the  plane  of  the  slab- 

The  parameters  S  and  5  denote,  respectively,  the  angle  between 
n  and  ex  and  %  and  e . .  In  terms  of  the  unit  vectors  cx,  ey,  and 
e_,  H  and  f  have  the  following  representations: 

n  =  cos  8  e  -  sin  8  e  (2-1) 

x  v 


4 


sin  0  o  +  cos  0  e  )  sin  6  +  cos  6  e  . 
x  y  z 


(2.2) 


The  unit  normal  pointing  into  the  illuminated  surface  of  the  slab 
will  be  denoted  by  n  and  is  given  by: 


n  =  £  x  n 


n  =  cos  6  (sin  6  ev  +  cos  0  e  )  -  sin  6  e  . 

A  >  £ 


(2.3) 


The  localizer  antenna  has  coordinates  (0,0  H)  relative  to  the  x, 
y,  z  coordinate  system,  the  height,  H,  o_  the  antenna  above  the 
ground  is  typically  about  10  feet.  The  midpoint  of  the  base  of 
the  slab  has  coordinates  (x^.y^.z^),  z^  being  the  distance  above 
the  ground  plane  of  the  base  of  the  slab.  The  angle  1 1>  is  just  the 
azimuth  angle  of  the  midpoint  (x^y^Zj): 


<i>  =  tan"1  . 


(2.4) 


Finally,  the  point  with  coordinates  (x2,y2,z2)  represents  the 
location  of  the  localizer  receiver. 

In  the  physical  optics  approximation,  the  scattered  electric 

field  £  in  the  Fraunhofer  or  radiation  zone  of  the  slab  is  given 
s 


where 


=  —  M 

2 it  \e  ) 


1/2  elkR2 


R2x^R2  x  ,  (2.5) 


J  =  /  (nxHj)  e"ik^2^)ds  . 


(2.6) 


In  F.quations  (2.5)  and  (2.6),  k  =  2ti/X  is  the  wavenumber  of  the 
radiation  (A  is  the  wavelength),  and  p  and  e  denote,  respectively, 
the  permeability  and  permittivity  of  space.  The  unit  vector  R2 
points  from  the  midpoint  of  the  base  of  the  slab  toward  the 
receiver  point: 


R2  =  ^2^R2  ’ 


(2.7) 
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where 


R2  n  (x2'xl)  ex  +  (y2"ylJ  ey  +  -22*  21^  ez 


(2.8) 


and 


R2  =  [(-VX1)2  +  ^2*yl)2  +  (Z2'Z1)2] 


1/2 


(2.9) 


The  integral  J  defined  by  Equation  (2.6)  is  carried  out  over  the 
illuminated  surface  of  the  slab.  The  vector  n  is  the  unit  normal 

■f 

pointing  into  the  surface  of  the  scatterer  while  H-  is  the  incident 

•y  1 

magnetic  field.  The  vector  r  is  the  position  vectoi  of  an  arbi¬ 
trary  point  on  the  illuminated  surface  relative  to  the  point  (x^, 
y^,Zj).  A  detailed  derivation  of  Equations  (2.5)  and  (2.6)  can  be 
found  in  Section  2  of  Reference  1. 

It  was  shown  in  Reference  1  that  the  incident  magnetic  field 
IL  at  a  field  point  (x,y,z)  produced  by  a  horizontally  polarized 
localizer  antenna  and  its  image  in  the  ground  plane  can  be  repre¬ 
sented  approximately  as  follows: 


Hi  = 


;2  0) 


1/2 


ikR 


Eo  V 


1  -  e 


2ikzH/Dp 


(2.10) 


In  Equation  (2.10),  Eq  is  an  amplitude  associated  with  the  locali¬ 
zer  antenna  while  f (4>)  is  the  horizontal  field  pattern  of  the 
antenna  (<*>  =  tan  *  (y/x))  .  The  parameters  R  and  Dp  denote, 
respectively,  the  slant  distance  and  horizontal  distance  from  the 
localizer  to  the  field  point: 


R  -  [x2  +  y2  +  (z-H)2]17"  , 

r  2  2I ^2 

Dp  -  [x2  .  y2J 


(2.11) 

(2.12) 


Equation  (2.10)  provides  an  accurate  representation  of  Hj  at 
field  points  satisfying  the  so-called  small  elevation  angle  condi¬ 
tion: 


:/Dp  «  1 


(2.13) 


It  will  be  assumed  in  the  present  work  that  the  inequality  (2.13) 
is  satisfied  by  all  points  on  the  surface  of  the  slab. 

The  vector  R  drawn  from  the  localizer  antenna  to  a  point  P 
on  the  surface  of  the  slab  can  be  represented  as  follows: 


■+  ->• 
R  =  R1  +  r 


(2.14) 


where  R^  =  ex  +  y^  ey  +  (z^-H)  e?  is  the  position  vector  of  the 
the  midpoint  of  the  base  of  the  slab  relative  to  the  antenna  and 
r  connects  the  midpoint  of  the  base  to  the  point  P.  To  first  order 
in  r/R^,  R  =  | R I  can  be  represented  as  follows: 


R  =  R^  +  ft.-r 


(2.15) 


where  R^  =  R^/R^.  The  vector  r  can  be  expressed  in  terms  of  the 
unit  vectors  n  and  £  as  follows: 


r  =  n  h  +  £  £  , 


(2 • 16 J 


where  n  ranges  from  -  L/2  to  1/2  (L  being  the  length  of  the  base 

of  the  slab)  and  £  ranges  from  0  to  h  (h  being  the  slant  height  of 

the  slab).  Using  Equations  (2.1)  and  (2.2)  to  express  n and  £  in 

terms  of  ev,  e,,  and  e„,  Equation  (2.16)  can  be  written  as 
x  >  *• 

follows : 


r  =  (n  cos  8  +  £  sin  6  sin  6)  e  +  (£  cos  6  sin  &  -  n  sin  0)e 

x  y 

+  £  cos  o  e  .  (2-17) 

*■'  A 

For  small  angles  of  elevation,  the  unit  vector  Rj  can  be  written 
approximately  as 


R.  =  cos  ip  e  +  sin  ip  e  + 


Uj-H)  . 


(2.18) 


where  Dp^  =  ^Xj  +  yj  and  is  the  azimuth  angle  of  the  midpoint 
(cf.  Equation  (2.4)).  Substituting  Equations  (2.17)  and  (2.18) 
into  Equation  (2.15)  yields  the  following  approximate  expression 
for  R: 
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iM'n."  ■  , . .  . . in 


R  ;  Rj  +  n  cos(6+$)  +  £  sin  6  sin 


cosS  (z.-H) 

(0+«)  +  € - n—1 -  (2-1! 


Introducing  the  angle  of  incidence  y  =  */2  ~  (9+v0  *  Equation 
(2.19)  can  be  rewritten  more  compactly  as  follows: 


l  «  Rj  <■  n  sin  v  *  ?  s 


coso  (z.-H) 
in  6  cos  y  +  C  - n - 


(2.20) 


The  angle  y  is  just  the  angle  between  the  projection  of  Rj  onto 
the  x-v  plane  and  the  projection  of  the  unit  normal  n  onto  the  x-y 
plane . 

The  approximate  expression  for  the  distance  R  given  in  Equa¬ 
tion  (2.20)  can  now  be  substituted  into  the  appropriate  complex 
exponential  in  Equation  (2.10).  Some  additional  approximations 

•t 

will  now  be  made  to  render  the  expression  for  H.  gi'sn  in  Equation 
(2.10)  more  suitable  for  the  integration  indicated  in  Equation 
(2.6).  Specifically,  ue  will  assume  that  is  so  much  greater 
than  L,  the  length  of  the  slab,  and  that  f(£)  is  sufficiently 
slowly  varying  then  f($)  can  be  replaced  by  f(i^)(£  -  tan’ 1  (y  /x } )) 
for  all  points  on  the  surface  of  the  slab.  That  is  to  say,  we  will 
assume  that  the  horizontal  angle  subtended  by  the  slab  at  the  lo¬ 
calizer  antenna  is  so  small  that  the  variations  in  f(4>)  over  the 
surface  of  the  slab  can  be  ignored. 

In  like  manner,  it  will  be  assumed  that,  to  a  good  appro¬ 
ximation,  the  amplitude  factor  R~*  in  (2.10)  can  be  replaced  by 
Rj1  and  that  the  parameter  in  the  expression  in  parentheses  can 
be  replaced  by  Dpj*  Incorporating  all  of  these  approximations 
into  Equation  (2.10),  we  obtain  the  following  expression  for  H. 
on  the  surface  of  the  slab:  . 


'st1'!  A  hi! 


5.  =  2ie  (*)  E  f(*)eikRl  eikn  sin  Y  •  e1 
1  zW  0  R 


k  ?sin  6  cos  y 


iks,  [H+C  cos  S]/D  . 

>  1  pl  •  sin  g—- 1 z j +  C  cos  6]  . 


(2.21) 


Consider  now  the  vector  R7  from  the  point  1  * “l^  t0  the 

observation  point  »2?)  •  h’e  will  denote  by  D^-,  the  horizontal 

distance  between  these  two  points: 


.2  '  [CVX1>2  *  ‘>-2->’l)2]1/2 


Assuming  that  the  elevation  angle  of  the  vector  R7  is  very  small 

A  “ 

(i.e.  Dj)7  >>  Js?-ZjJ),  the  unit  vector  R2  =  R7/R7  takes  on  the 
approximate  form: 


R_  =  cos  c  -  sin  o  e..  +  — £ — —  e_  , 
2  x  y  D  4 


(22‘2i^  - 


P2 


(2.22) 


where  0  is  the  angle  between  the  projection  of  R7  onto  the  x-y 
plane  and  the  x-axis.  Consequently,  the  quantity  ft7-r  appearing 
in  Equation  (2.6)  can  be  written  approximately  as  follows: 

„  _  5cos  6  (Z7-Z.) 

R7-r  =  n  cos  (0-6)  -  5  sin  6  sin  (o-6)  +  - ~ — “ - ’  (2.23) 

V 

where  use  has  been  made  of  Equation  (2.17).  Introducing  the  aspect 
angle  g  =  «/ 2  -  (0-6),  Equation  (2.25)  can  be  rewritten  more 
compactly  as: 


R2*r  =  n  sin  fc  -  K  sin  6  cos  S  + 


C  cos  5  (z7-z.) 


(2.24) 


The  angle  S  is  just  the  angle  between  the  projection  of  R7  ont^ 
the  x-y  plane  and  the  projection  of  (-n)  onto  the  x-y  plane  where 
n  is  the  unit  inward  normal  to  the  surface. 
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Th?  int~i;ral  in  Equation  (2.6)  can  now  be  evaluated  quite 

A  *t  A  ^ 

readily.  Substituting  the  expressions  for  n,  rK  and  R2*r  8iven 
respectively,  in  Equations  (2.3)  and  (2.21),  and  (2.24)  and  inte¬ 
grating  with  respect  to  n  (-L/2  <  n  <  +  L/2)  and  C  (05€5h) ,  we 
obtain  the  following  result  for  the  scattered  field  Eg: 

ik(R1+R2) 

E  =  -  0Qf C'P)  J["R~  cos  5  sinc  (sin  y  *  sin 

X  2 


[Sjxtftjxn)] 


(2.25) 


where  sinc(x)  =  sin  x/x  and  Ij  and  I2  are  defined  as  follows: 


eiAh  -  1 


02ikHz-i/Dnl  [iBh  -1 
t  _  e  i  pi  Le  -  jj 

lZ  B 


(2.26) 


(2.27) 


where 


A  =  kjsin  6  (cos  6  +  cos  6)  +  cos  6 


B  =  A  ♦  ~  cos  6 
Pi 


(2.29) 


For  the  small  angles  of  elevation  considered  here,  Eg  is  essentially 
parallel  to  the  ground  and  its  magnitude  Es,  is  given  by: 


ik(R1+R2) 

Es  =  ‘  ll  Eof^  tj-Rj 


kL 

cos  0  cos  8  sine  y-(sin  y  -  sin  8) * 


■M  • 


(2.30) 


To  take  into  account  the  image  of  the  slab  in  the  ground 
plane,  the  expression  for  Eg  given  in  Equation  (2.30)  should  be 
evaluated  at  the  image  of  the  observation  point  (x2,y2,-s2)  and 


the  result  subtracted  from  Equation  (2.50).  The  final  result  for 
the  total  scattered  field  E,.T  (direct  plus  ground  reflected!  at  the 

w  I 

receiver  is  just: 


^ 71^7)  ”  E^ (x^ ^ ,  * .. 7)  , 


(2.51) 


where  is  given  in  Equation  (2.30). 

The  expressions  for  the  scattered  electric  field  generated 
by  an  elevated,  tilted,  rectangular  slab  developed  in  this  section 
constitute  the  basic  mathematical  formulas  used  in  our  localizer 
computer  program  for  treating  ILS  signal  scattering  by  structures 
with  plane  surfaces. 


2.3  LOCALIZER  SIGNAL  SCATTERING  BY  A  RIGHT,  CIRCULAR  CYLINDER 

In  this  section,  localizer  signal  scattering  by  a  right, 
circular  cylinder  is  investigated.  This  problem  is  of  practical 
importance  since  cylindrical  structures,  such  as  gas  storage 
tanks,  are  frequently  to  be  found  near  airport  runways. 

The  geometry  of  the  problem  is  illustrated  in  Figure  2.2. 

The  x-y  plane  is  the  ground  plane  which  we  assume  to  be  perfectly 
flat  and  perfectly  conducting.  The  x-axis  represents  the  center 
line  of  the  runway  while  the  z-axis  is  perpendicular  to  the  ground 
plane  and  points  out  of  the  page.  The  unit  vectors  in  the  x,  y, 
and  z  direction  are  denoted  by  ev,  ey,  and  6_,  respectively.  The 
coordinates  of  the  localizer  relative  to  the  x,y,z  coordinate 
system  are  (0,0,H).  The  cylinder  is  assumed  to  be  perfectly  con¬ 
ducting  and  to  be  resting  on  the  ground.  The  radius  of  the  cylin¬ 
der  is  denoted  by  a.  The  coordinates  of  the  center  of  the  base 
of  the  cylinder  are  denoted  !<y  (Xj,Vj,0)-  The  azimuth  angle  of 
this  center  point  is  denoted  by  £(«=tan~ 1  (y^/Xj)  )  .  The  point  with 
coordinates  (x,,y7,z7)  represents  the  location  of  the  localizer 
receiver.  The  vectors  R^  -*nd  R?  denote,  respectively,  the  vector 
connecting  the  localizer  enna  to  the  point  (Xj,y^,0)  and  the 
vector  connecting  the  point  (x^y  0)  to  the  receiver  point 

(x?  ,y7,i_7). 


"W'lfiMi  i  * 


R1  =  X1  ex  +  yi  ev  '  H  e: 


R,  =  (x2-x1)  ex  +  (Vj-Vjj  ey 


(2.32) 

(2.53) 


Finally,  the  angle  3  is  just  the  obtuse  angle  between  the  projec- 
tions  of  Rj  and  R?  onto  the  x-y  plane. 

Equations  (2-5)  and  (2.6)  will  be  used  to  calculate  the 
scattered  electric  field  in  the  radiation  tone  of  the  cylinder. 

The  vector  r  appearing  in  Equation  (2.6)  should  now  be  interpreted 
as  connecting  the  center  point  (Xj,y^,0)  to  an  arbitrary  point  P 
with  coordinates  (x,y,t)  eft  the  illuminated  surface  of  the  cylinder. 
Clearly,  the  vector  r  can  be  represented  as  follows: 

r  =  a  cos  n  +  a  sin  4  Z  +  z  e_  ,  (2.34) 


where 


n  =  cos  f  ex  +  sin  f  e^_  ,  (2.35) 

and 

t~~  sir.  #  e  +  cos  $  e  .  (2.36) 

The  vector  n  is  the  unit  vector  in  the  direction  of  the  projection 
of  Rj  onto  the  x-y  plane  and  the  angle  4  is  angle  between  n  and 
the  projection  of  r  onto  the  x-y  plane.  Since  we  are  assuming  that 
the  cylinder  is  in  the  far  field  of  the  localizer,  o  ranges,  to 
good  approximation,  from  b/2  to  Ss/Z  on  the  illuminated  surface 
of  the  cylinder.  Substituting  Equations  (2.35)  and  (2.56)  into 
(2.34),  we  obtain  the  following  expression  for  the  position  vector 
r: 


r  =  a  cos(^+f)  e  +  a  sin(£+f)  e  +  z  e_  .  (2.57) 

A  )  ^ 

As  in  Section  2-2,  we  will  make  the  small  elevation  angle 
approximation  so  that,  to  a  good  approximation,  the  total  incident 
magnetic  field  H-  (direct  plus  ground  reflected)  at  a  point  P  on 
the  surface  of  the  cylinder  is  given  by: 
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h.  =  -  e.l  — 


lkRl  f  2ikzH/D  ,1  ikR. . r 

E  ff*)  1  -  e  pi  e  1 

O  '  Kj 


P1  e  1 


(2.3S 


*  -r  0  7  \jO 

where  Rj  =  Rj/R^  and  ^pl  =  ‘x£  +  >'jl  *  For  snail  angles  of 
elevation,  R^  can  be  represented  approximately  as  follow*- 

z.  *  .  a  h 

R.  =  cos  e  +  sm  ^  e  -  jt —  e  (2.3P 

i  x  >  u{)1  ^ 

Consequently,  the  inner  product  R^*r  becomes  simply  (cf.  Eq. 
(2.57)3: 


Rrr  =  a  cos * -  f 


(2.40 


Substituting  Equation  (2.40)  into  (2.38),  we  obtain  the  following 


result  for  H. 


-  V\1/2 

Hi =  2i  Mf) 


KofW  57  el 


ika  cos$  .  kHz 


dpi  ' 


(2.41 


The  unit  inward  normal  to  the  cylindrical  surface  at  the 
point  P  is  given  by: 

n  =  -  cos  (6+9)  e  -  sin  ($+$)  e  .  (2.42 

x  ; 

A  ■+ 

For  small  angles  of  elevation,  the  unit  vector  R7  =  R7/R?  can 
be  written  approximately  as  follows: 


R,  ;  cos  v  e  -  sin  y  e  +  jr=—  e  , 

x  -  Up2  i 


(2.45 


where  =  ^(x-j-Xj)"  +  (y7-y^)“J  ~  and  y  is  the  angle  between  th 
x-axis  and  the  projection  of  R,  onto  the  x-y  plane.  Consequently 
the  inner  product  R?*r  can  be  expres'.-ed  as  follows: 


R->-r  =  a  cos  (0+$+y3  * 


(2.44 


Figure  2-1.  Geometry  of  the  Cylinder  Scattering  Problem 


Referring  to  Figure  2-2,  we  note  that  the  following  identity 
holds  among  the  angles  #,  S,  and  y: 

Cf*f)  =  S  -  B  .  (2-4 

Substituting  Equation  (2-4S)  into  (2.44),  we  obtain  the  following 
results: 


R?-r  =  --a  cos  (£-S) 


(2-4 


The  scattering  integral  (2.6)  can  now  be  evaluated  quite 

•W  A  A  ^ 

readily  by  using  the  expressions  for  Hj,  n,  and  R?-r  given,  re¬ 
spectively.  in  Equations  (2.41),  (2.42),  and  (2.46).  The  inte¬ 
grations  are  performed  with  respect  to  #  (s/2  <  ♦  <  3v/2)  and 
z  (-h  <  z  <  +  h).  The  parameter  h  is  the  height  of  the  cylinder 
and  the  range  of  the  z  integration  is  fro*  -h  to  +h  to  take  into 
account  the  image  of  the  cylinder  in  the  ground  plane.  We  find 


that  for  sc3ll  angles  of  elevation  (D^  >5,  -?)  ,  Eg  is  polarized 
essentially  parallel  to  the  ground  plane  and  that  Its  nagnitude, 
E_,  Is  given  by: 


ikiR,-R,) 


Eg  =  E Qf  (?)  jr^ 


a  h  sinc  1  hfcr ' 
- 


f  H  2,  1 

sint  k  hj  gp—  +  k3—  ‘  I  (S)  » 

L  ?1 


(2.47) 


vnere 


7'"’  ~  2Ika  cos  3/2  cos  C#~S/2) 

cos  (e-£)  e  d#  .  (2.4S) 

-/? 

The  integral  1(3)  can  be  expressed  as  an  infinite  scries  of  Bessel 
functions  as  follows: 

ICS)  =  -  -  i  cos  3/2  Jj(v)  +  2  cos”  3/2  JQ(v} 


sin  3  ...  , , 

S3  (l+COSS)  S1S  {‘v3  S1R  sj 


n,->  V  7  COS  f2n-l)S/2 

cos  s/2  ^  J2n(v)  - 2^1 - 


cos  { 2n+l)6/2l 

2r.+l  J  1 

(2.49) 


shore  v  =  2ka  cos  3/2  and  J  is  the  Bessel  function  of  order  n. 

n 

A  very  efficient  subroutine  has  been  developed  for  calculating 
1(3)  to  any  desired  degree  of  accuracy  using  the  above  series 
representation. 

It  should  be  noted  that  Equation  (2.47)  represents  the  physical 
optics  approxisation  to  the  scattered  field  produced  by  a  right, 
circular  cylinder  and,  consequently,  is  only  valid  when  a  »  A 
(ka  »  1), 
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To  illustrate  the  predictions  of  our  cylinder  model,  we  have 
plotted  in  Figures  2-3  and  2.4  the  microampere  deviations  (static 
and  dynamic)  along  the  centerline  of  a  runway  caused  by  a  cylinder 
80  f  et  in  diameter  and  50  feet  high  whose  center  is  located  5260 
feet  down  the  runway  from  a  V-Ring  localizer  and  470  feet  off  the 
runway  centerline.  The  abscissa  in  each  figure  represents  dis¬ 
tance  from  the  localizer.  The  aircraft  is  assumed  to  be  executing 
level  flight  at  an  altitude  of  50  feet. 
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3.  GLIDE  SLOPE  TERRAIN  SCATTERING 


3-1  INTRODUCTION 

The  majority  of  ILS  glide  slope  antenna  arrays  currently  in 
use  are  o£  the  image  variety.  Consequently,  the  quality  of  the 
guidance  information  provided  by  these  arrays  is  strongly  dependent 
upon  the  characteristics  of  the  surrounding  terrain.  In  this 
report,  a  mathematical  model  is  developed  for  predicting  the  dis¬ 
tortions  in  glide  slope  antenna  patterns  produced  by  certain  types 
of  large  scale  terrain  irregularities.  Results  are  presented 
illustrating  glide  slope  performance  for  some  idealized  terrain 
configurations . 

5 . 2  THEORY 

The  vertical  antenna  pattern  of  an  image-type  glido  slope 
array  is  determined  by  the  interference  between  the  direct  and 
ground  reflecteu  radiation.  The  ground  reflected  or  scattered 
radiation  is  generated  by  the  charges  and  currents  induced  in  the 
ground  by  the  incident  radiation  from  the  glide  slope  array.  For 
simplicity,  we  will  assume  that  the  ground  plane  is  perfectly 
conducting.  For  a  perfect  conductor,  the  surface  current  density 
K  is  given  by: 


K  =  n  x  H 


(3.1) 


where  n  is  the  unit  normal  pointing  out  of  the  ground  and  H  is  the 
cotal  magnetic  field: 


H  =  H.  +  H 
x  s 


(3.2) 


where  H.  denotes  the  incident  and  H  denotes  the  scattered  magnetic 
field. 

Given  the  surface  current  density  K  on  the  ground,  the  scat- 
-> 

tered  field  H.  at  a  receiver  above  ground  is  given  by  the  following 
surface  integral: 
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r  ■1«I,|  !'1J,|I(|H'V  !,ll  I1  1  lf|^'^^nirpii(WiMi[#T)l1f[l,irt||fiilHi|  inn  "Hit 


(3.3) 


■> 

H_  (r 


S-hf* 


K(r)  x  VG  (r1,r)  ds 


Here  and  r  denote,  respectively,  the  position  vector  of  the 
receiver  and  the  position  vector  of  an  arbitrary  source  point  on 
the  ground  plane  S  relative  to  some  origin  of  coordinates.  The  two- 
point  Green’s  Function  G(r^,r)  is  given  by 


G  (tj.r)  = 


ik|rrr[ 


(3.4) 


where  k  =  2ts/\  and  X  is  the  wavelength  of  the  incident  radiation. 

To  specify  the  current  density  K,  we  will  adopt  here  the  same 
physical  optics  model  used  in  Reference  1  for  treating  high  fre¬ 
quency  scattering  by  perfect  conductors.  Specifically,  we  will 
assume  that  on  those  portions  of  the  ground  not  directly  illuminat- 
ed  by  the  array,  K  is  identically  zero  and  that  on  those  portions 
of  the  ground  which  are  directly  illuminated,  K  is  given  by  twice 
the  tangential  component  of  the  incident  magnetic  field: 


K  =  0  on  S_  , 

K  =  ZfnxHp  on  S+  , 


(3.5) 


(3.6) 


where  S_  and  S+  denote,  respectively,  the  unilluminated  and  illu¬ 
minated  portions  of  the  ground  plane  S.  The  two  types  of  regions 
are  illustrated  in  Figure  3-1.  The  current  distribution  given  in 
Equations  (3.5)  and  (3.6)  should  be  reasonably  accurate  when  the 
characteristic  dimensions  of  the  terrain  irregularities  are  large 
compared  to  the  wavelength  X.  Substituting  Equations  (5.5)  and 
(3.6)  into  Equation  (3.3),  we  obtain  the  following  approximate 
expression  for  the  scattered  magnetic  field  at  the  receiver  point 


(rp  »  f  [n(r)  x  H± (r)J  x  ^G  (r^r 


)  ds  . 


(3.7) 


21 


Antenna 


Figure  3-1.  Cross  Section  of  Terrain  Showing  Illuminated 
and  Unilluminated  Regions 

In  general,  the  evaluation  of  the  integral  in  Equation  (3.7) 
for  an  arbitrarily  irregular  ground  plane  is  a  prohibitively  time- 
consuming  task  even  for  the  most  rapid  digital  computer.  However, 
terrains  characterized  by  irregularities  which  vary  only  along  one 
coordinate  axis  can  be  treated  with  relative  ease.  Consequently, 
we  will  confine  our  attention  at  present  to  such  one-dimensional 
ground  planes.  Specifically,  the  touchdown  point  on  the  center- 
line  of  the  runway  will  bt  chosen  as  the  origin  of  coordinates  0. 
The  z-axis  is  chosen  to  be  the  vertical  axis  passing  through  0  and 
pointing  out  of  the  ground  while  the  x-  and  y-axes  lie  in  the 
horizontal  plane  whose  normal  parallels  the  z-axis.  The  x-axis 
points  along  the  centerline  of  the  runway  while  the  y-axis  is 
perpendicular  to  the  centerline.  The  one-dimensional  irregular¬ 
ities  are  assumed  to  vary  only  with  x.  That  is  to  say,  the  equa¬ 
tion  of  the  surface  of  the  ground  is  of  the  form  z  =  f(x)  and  is 
independent  of  variations  in  y. 

Consider  now  the  problem  of  a  horizontal  dipole  mounted  above 
such  a  ground  plane  and  oriented  parallel  to  the  y-axis,  the 
invariant  axis  of  the  terrain  irregularities.  This  problem  is  of 


fundamental  importance  since  the  majority  of  glide  slope  antenna 
systems  consist  of  arrays  of  such  dipoles.  The  geometry  of  a 
typical  problem  is  illustrated  in  Figure  3-2.  In  Figure  3-2,  the 
y-axis  points  into  the  page.  Since  the  terrain  is  assumed  to  be 
y-invariant,  the  entire  griund  plane  can  be  generated  by  project¬ 
ing  the  cross  section  depicted  in  Figure  3.2  in  the  y-direction 
from  -®  to  +<*>. 

•f 

To  calculate  the  scattered  field  Hg,  generated  by  a  given 
terrain  contour,  the  ground  surface  is  first  divided  into  a  number 

of  connected  planar  sections  which  extend  from  -»  to  +®  in  the 

2 

y-direction.  The  side  view  of  such  a  section  >s  represented  by 
the  darkened  segment  of  the  contour  depicted  in  Figure  3-2.  The 

points  (xoi»2oi^  anc^x02’z02^  rePresent  the  end  points  of  the 
section.  The  length  L  of  a  given  segment  is  so  chosen  that  the 
unit  normal  n  is  essentially  constant  for  x^j  <  x  <  x^  and 


Figure  3-2.  Geometry  of  a  Typical  One-Dimensional  Ground 
Reflection  Problem 
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L  ~  [*-x02~X01^  +  ^202'201^~] 


?1l/2 


(3.8) 


Once  the  contour  has  been  segmented  into  chese  planar  sections. 
Equation  (3.7)  can  be  applied  to  each  section  and  the  total  scat- 
tered  field  Hg  can  then  be  obtained  by  summing  the  contributions 
from  all  the  sections.  Ke  will  therefore  concentrate  upon  obtain¬ 
ing  an  expression  for  the  contribution  to  the  total  scattered 
field  Hs  from  a  typical  planar  segment  such  as  the  one  depicted 
in  Figure  3.2.  It  will  be  assumed  that  only  those  segments  which 
lie  between  the  transmitter  and  receiver  contribute  to  H  . 

C 

Let  hg  denote  the  contribution  to  Hg  from  the  planai  section 
shown  in  Figure  3.2: 


hs (V  =  Jr  J [n  x  H± (r)J  x  ?G  (rj.r)  ds  , 


(3-9) 


where  the  integral  is  carried  out  over  the  surface  of  the  segment 
(cf.  Eq.  (3.7)).  Let  e^,  ev,  and  e„  denote,  respectively,  the 
unit  vectors  in  the  x,  y,  and  z  directions.  In  terms  cf  these 
unit  vectors,  the  unit  normal  n  is  given  by: 


-  sm  e  e„  +  cos  c  e 

X  4 


(3.10) 


where  e  is  the  angle  between  the  planar  section  and  the  x-axis 

(cf.  Fig.  5.2).  The  x,z  coordinates  of  the  transmitting  dipole 

and  the  receiver  are  denoted  by  (0,h)  and  (x.,tj),  respectively. 

For  simplicity,  we  will  assume  that  the  receiver  is  above  the 

centerline  of  the  runway  so  that 

y ^  =  0 .  (3-11) 

The  V  displacement  of  the  center  of  the  transmitting  dipole 
will  be  denoted  by  y  .  Glide  slope  arrays  are  generally  located 
several  hundred  feet  from  the  centerline  of  the  runway. 

The  incident  magnetic  field  H.  produced  by  the  transmitting 
dipole  is  given  in  terms  of  the  volume  current  density  J  in  the 
dipole  by  the  following  formula  (Reference,*  Equation  (2.23)). 


H^r) 


=  L.  J(d(r’)  x  V'Gtf.r’^dv*  , 


(3.12) 


where  V'  operates  on  the  primed  coordinates.  For  simplicity,  we 
will  assume  that  all  dimensions  of  the  dipole  are  small  compared 
with  the  wavelength  X.  This  assumption  is  not  strictly  valid  since 
the  dipoles  which  make  up  most  glide  slope  arrays  are  X/2  in 
length.  However,  the  assumption  of  electrical  smallness  simplifies 
the  calculations  enormously  and  besides,  the  exact  form  of  the 
dipole's  azimuthal  radiation  pattern  is  not  really  relevant  since 
we  are  concerned  here  with  the  vertical  radiation  patterns  of 
glide  slope  arrays.  Consequently,  we  will  assume,  for  simplicity, 
that  the  transmitter  is  a  point  dipole: 


J(r')  =  Jo  s  o(x')  <5(y'-ya)  o(z'-h) 


(3.133 


where  J  is  some  constant  and  as  noted  previsously,  v  is  the 
o  a 

displacement  of  the  dipole  from  the  centerline  of  the  runway.  The 
S- functions  are  just  the  Dirac  delta  functions.  We  will  make  the 
further  assumption  that  the  planar  section  is  many  wavelengths 
distant  from  the  transmitter  so  that  to  a  good  approximation, 

*t  **  + 

V’G  (r,r*)  can  be  written  as  follows: 


V'G  (r ,r' )  « 


.k  (r-r'J  gifcjr-r* 
ir-f*!  !r-r’ I 


(5.14) 


Substituting  Equations  (3.13)  and  (5.14)  into  Equation  (5.12),  we 
obtain  the  following  expression  for  the  incident  magnetic  field 

iyr): 

-  ikJe  r  „  -j 

Hi<r>  *  TT  V"  [X  6z  '  (Z'h)  C*J  ’  C5-15J 

where  r  =  x  e  +  y  ev  +  z  c_  and  D  is  defined  as  follows: 


=  [x2  *  (y-ya)2  +  (s-h)2]' 


(3.16) 


Equations  (5.10)  and  (5.16)  can  now  be  substituted  into 
Equation  (3.9)  and  the  magnetic  field  hg  at  the  receiving  point 
Tj  =  Xj  cx  +  Zj  e_  due  to  the  single  planar  section  of  terrain 
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evaluated.  It  will  be  assumed  that  the  receiver  is  many  wavelengths 
distant  from  the  terrain  segment  so  that  ?G  (r^,r)  can  be  written 
approximately  as  follows: 


7C-  (r. ,r)  =  -  ik  ■  ‘  - 


(?,-?)  „ik|V?i 


|rrr|  IV7! 


(3.1?) 


Actually,  only  the  z-cosponent  of  hg  need  be  calculated  since 
glide  slope  receivers  respond  primarily  to  this  component.  The 
z-  component  of  h,  will  be  denoted  by  hs_.  Substituting  Equations 
(3.10),  (3.15),  and  (3.17)  into  Equation  (5.9),  we  obtain  the 
following  expression  for  h  at  the  receiver: 


t  ,■*  ,  r2jo  f  Cx1-x3  ex  Sin  c  -  (z-h)  cos  e)  ik(r#+R)j 

h  (r.)  - y-  I - ? — 5 -  e  k  'ds,  (o.l8) 

Si  1  8-“  J  C"  R~ 

where  R  =  ^(Xj-x)'  +  >'”  +  (Zj-z)^p^.  The  element  of  area  ds  in 
(5.18)  is  equal  to  dv  dq,  where  n  is  a  variable  ranging  from  0  to 
L  along  the  surface  of  the  planar  section  of  terrain  perpendicular 
to  the  y-axis-  The  x~  and  y-coordinates  of  a  point  on  the  surface 
of  the  planar  section  can  be  expressed  in  terms  of  n  as  follows 
(cf.  Fig.  3-2): 


x  =  Xq j  +  n  cos  e  ,  0  <  n  <  L  , 


q  sin  z  ,  0  <  n  <  L  , 


(5.19) 


(3.20) 


Substituting  Equations  (3.19)  and  (3.20)  into  (3.18),  we  find  that: 

k  r  I 

hsz(rl^  =  — 2~  L(20l"h)  cos  c  '  *01  sin  £J 


f  f  {xrx0l'  ncos  £)  ik(D+R) ,  ,  , 


where  D  and  R  are  expressed  as  functions  of  the  variables  of  inte¬ 
gration  y  and  n  as  follows: 


=  [(x01+n  cos  0)2  +  (>'ya)2+  (20i'h+n  sir*  e-2]  1  •  C3*22) 


R 


r  9  2  2i 1/2 

=  ^Xj-XQj-n  cos  £)"  +  y  +(z1-z01-n  sin  e)  .  (3.25) 

The  integration  with  respect  to  y  can  be  performed  approximate¬ 
ly  using  the  method  cf  stationary  phase  leaving  only  an  integral 

■y 

with  respect  to  n  to  be  evaluated  numerically :~ 


h  (r.)  = 

C  T  V  1  * 


k2J 


*sil*l'  "  “ T 

ft  u 


[<*orh) 


cos  e  -  sm  e 


1  >1 / 2  i* 

J  A  ° 


r'/  4 


I  ,  (3-24) 


where 


I  =  f  eikFCll)  G  (n)  dn  , 


(3.23) 


and  F  and  G  are  defined  as  follows: 


A  = 


B  = 


C  = 


[(xrxorn  cos  E)2  +  Czrz0rn  sin  £)j 

r  9  ,-ji/2 

[(x0I+n  cos  e)“  *  (zQ1-h+n  sin  e)“J 

v2  11/2 


,,1/2 


1  ♦ 


'  a 


(A+B)  ‘ 


,  D  =  CB,  R  =  CA, 


■■?(«*•) 


F  =  C  (A+B)  , 


Xj  -  x0j  -  n  cos  c 
02  R2  S1'2 


(3-26) 


The  z-component  of  the  total  scattered  magnetic  field  at  the 


receiver  is  obtained  by  summing  up  the  contributions  h  from  all 


of  the  illuminated  planar  sections  lying  between  the  transmitter 
and  receiver.  Symbolically,  this  result  can  be  expressed  as 
follows : 


H 


sz 


-I 


(3.27) 


The  contribution  h  from  each  section  is  calculated  using  Equa¬ 


tions  (3.24),  (5.25)  and  (3.26),  the  integral  in  (3.25)  being 
performed  numerically. 

The  z-component  of  the  total  magnetic  field  (direct  plus 
scattered)  at  the  receiver  is  given  by  (cf.  Equation  (5.15)): 


ikJ 


ikR, 


H 


o  ^  v  a  u 

71  X1  +  Hs; 

R1 


(3.28) 


r7  7  7l  1  /  7 

x~  +  ya  +  (zj-h)~J  .  The  first  term  in  Equation 
(3.28)  is  just  the  z-component  of  the  expression  for  H-  given  ir. 
Equation  (3.15)  evaluated  at  the  receiver  point  (Xj.OrZj).  To 
calculate  the  vertical  radiation  pattern,  it  is  convenient  tc  use 


as  an  origin  of  coordinates  the  point  (0,y  ,0),  which  is  the  point 

3 


on  the  ground  directly  below  the  antenna.  Let  Rq  denote  the  dis¬ 
tance  from  this  point  to  the  receiver: 


R0 


-[4 


-  .?lm 


(3.29) 


We  will  assume  that  the  receiver  is  sufficiently  far  away  from  the 

-7 

transmitter  so  that  the  amplitude  factor  R."  in  (3.23)  can  be 

.  7  1 


replaced  by  R 


0 


Under  this  assumption.  Equation  (3.28)  can  be 


rewritten  as  follows: 


kJ 


ikR. 


H_  * 


o  e 


ikCR.-R,.) 


4v  R0  *lkR0 
+  | — f —  C 

kJo  X1 


H, 


1 

M 


(5.50) 
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The  expression  in  brackets  in  Equation  (3,30)  represents  the 
vertical  radiation  pattern  P(#)  of  the  dipole-terrain  system: 


PC#)  -  i  o 


ik(Kj-Rg)  4^  R0  *ikR0 

*  ET  it  e  «. 


(5.31) 


The  angle  #  is  just  the  elevation  angle  of  the  vector 
RQ  =  Xj  ex  +  y'a  ev  +  z j  e,  relative  to  the  x-y  plane: 


#  =  tan 


1  (v<*!  *  &'*} 


(3.32) 


For  an  array  of  such  dipoles,  one  would  calculate  a  P(#l  for  each 
elenent  and  then  sun  up  these  individual  patterns  tc  obtain  the 
composite  pattern  of  the  array-terrain  system. 


5.3  NUMERICAL  RESULTS 

In  this  section,  we  present  sose  of  the  preliminary  nunerical 
results  which  have  been  obtained  using  the  formal isa  developed  in 
the  preceding  section. 

In  Figure  3-5,  we  have  plotted  the  magnitude  of  the  sidehand- 
onlv  pattern  of  a  null  reference  glide  slope  array  as  a  function 
of  the  angle  of  elevation  d assuming  a  perfectly  flat  ground  plane. 
The  glide  slope  wavelength  A  was  set  at  3  feet  and  the  height  ft 
of  the  sideband  antenna  was  set  at  50  feet  (h=10A).  The  dashed 
curve  represents  the  ideal  sideband-only  pattern  based  upon  image 


tneory: 


(Pf  ~  2  | sin  (20-  sin  #-)  [ 


(3.35 


The  solid  curve  represents  the  sideband-only  pattern  calculated 
using  Equation  (5-31).  This  pattern  was  calculated  by  sussing  the 
contributions  from  the  first  5000  feet  of  ground  plane.  The  close 
agreement  between  the  two  curves  obviously  indicates  that  5000 
feet  of  flat  ground  plane  is  quite  sufficient  to  fora  the  null 
reference  SBO  pattern. 

As  a  further  Illustration  of  the  predictions  of  our  model,  we 
have  plotted  in  Figure  3-4  the  magnitude  of  the  sideband-only 
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igure  5-5.  Magnitude  of  Null  Reference  Sideband-Only 
Pattern 


----  iaage  Theory 


TSC  Model 


igure  5-4.  Magnitude  of  Capture  Effect  Sideband* 
Onsv  Pattern 
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pattern  of  a  capture- effect  array,  again  assuning  a  perfectly  flat 
ground  plane.  The  capture-effect  array  consists  of  three  dipoles 
placed  one  above  the  other  at  heights  which  were  set  at  15,  50, 
and  45  feet  (5X,  10X,  and  15X) .  All  three  dipoles  are  fed  sideband 
only  signal  with  the  following  relative  amplitudes:  (-.5,  1.0, 
-.5).  The  resulting  array  pattern  according  to  image  theory  is 
given  by: 


| P I  =  {2  sin(20r  sin  $)  [1  -  cos(10s  sin  9)11 


[5.34) 


The  dashed  curve  in  Figure  3.4  represents  the  ideal  inage  theory 
pattern  of  Equation  (5.34)  while  the  solid  line  represents  the 
pattern  obtained  by  applying  Equation  (3.31)  to  each  antenna 
element  and  then  summing  the  individual  patterns.  This  calculated 
pattern  represents  the  contributions  from  the  first  5000  feet  of 
ground  plane-  As  in  the  case  of  the  null  reference  SBO  pattern, 
the  agreement  between  the  two  curves  is  excellent. 

Finally,  to  illustrate  the  predictions  of  our  theoretical 
model  in  the  case  of  an  irregular  ground  plane,  we  applied  the 
model  equation  ef  Section  3.2  to  the  terrain  configuration  depicted 
in  Figure  3.5. 

transmitter 


Figure  3-5.  Sample  Siting  Proolcs 


Referring  to  Figure  3.5,  which  obviously  is  not  drawn  to  scale,  we 
note  that  the  terrain  in  front  of  the  antenna  consists  of  1200  feet 
of  flat  ground  terminating  in  a  40- foot  drop  to  a  lower  plateau 
which  then  extends  out  an  additional  3S00  feet.  This  type  of 
siting  problem  was  one  of  many  suggested  to  us  for  investigation 
by  the  FAA.  The  sideband  antenna  of  a  null  reference  glide  slope 
array  is  placed  at  a  height  of  30  feet  (10A)  above  the  first 
plateau  and,  therefore,  at  height  of  70  feet  (23. 3X)  above  the 
lower  plateau.  In  Figure  3-6,  we  have  plotted  the  null  reference 
S30  pattern  for  this  terrain  configuration  as  predicted  by  our 
model . 

The  pattern  depicted  in  Figure  3.6  exhibits  features  which 
reflect  the  two  effective  antenna  heights  (50  feet  and  70  feet) 
which  arise  because  of  the  40  foot  drop-off.  The  nulls  of  a  null 
reference  SBO  antenna  pattern  occur,  according  to  image  theory,  at 
those  values  of  the  elevation  angle  #  which  satisfy  the  following 
equation: 


sin  (kh  sin  6)  =  0 


(3.35) 


where  h  is  the  height  of  the  antenna  above  the  ground  plane.  For 
snail  values  of  6,  the  first  few  nulls  occur  at  the  following 
angles: 


HA 


,  n  «  0,  1,  2, 


(3.35a) 


For  X  =  3  feet  and  h  =  30  feet,  the  first  few  nulls  occur  at 
multiples  of  2.86’  while  for  h  =  70  feet  the  first  few  nulls 
occur  at  multiples  of  1.23°.  Referring  to  Figure  5.6,  we  note 
that  the  primary  ainisun  occurs  near  2.S6®,  the  glide  path  angle 
corresponding  to  the  50-foot  antenna  height.  In  addition,  there 
is  a  shallower  local  aininus  near  1.25s,  the  glide  path  angle 
corresponding  to  the  effective  70- foot  antenna  height.  Note  also 
that  there  is  an  inflection  point  near  2.46®  which  is  twice  the 
70- foot  glide  path  angle. 


In  suanary,  the  aain  features  of  the  pattern  depicted  in 
Figure  3.6  reflect  contributions  to  the  scattered  field  from  the 
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first  1200  feet  of  ground  plane  for  which  the  effective  antenna 
height  is  30  feet  (10a).  There  are,  however,  secondary  features 
which  clearly  reflect  contributions  to  the  scattered  field  froa 
the  3800  feet  of  ground  plane  beyond  the  drop-off  for  which  th  ■ 
effective  antenna  height  is  ?0  feet  (25. 5>.)  - 
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l\,  MODEL  OF  A  PARABOLIC  REFLECTOR  LOCALIZER  SYSTEM 


4.1  INTRODUCTION 

The  introduction  of  the  Texas  Instruments  (TI)  wide  aperture 
parabolic  reflector  localizer  raises  a  question  about  the  ability 
of  the  ILS1.0C  program-3  to  properly  predict  the  Course  Deviation 
Indication  (CDI)  for  arbitrary  signal  transmitting  systems.  The 
program  has  a  provision  for  accepting  measured  course  and  clearance 
signal  amplitudes  to  characterize  arbitrary  array  patterns.  This 
method  assumes  that  the  measured  fields  are  generated  by  linear 
arrays  elevated  to  a  discrete  height  above  the  ground  plane.  For 
such  a  case  there  should  be  no  anomalous  phase  variations  between 
the  carrier-plus-sidebands  and  sidebands-only  signal  components 
as  azimuth  ~ngle  is  varied. 

In  contrast,  the  parabolic  screen  of  the  TI  antenna  repre¬ 
sents  a  vertically-extended  source  with  a  vertical  pattern  which 
cannot  be  described  by  a  discrete  point  image  formula.  Moreover, 
it  is  not  clear  that  sidebands-only  and  carrier-plus-sidebands 
signals  are  generated  so  as  to  have  a  common  center  of  phase. 

These  considerations  raised  the  question  of  whether  the  use 
of  measured  TI  arrarv  patterns  without  detailed  phase  information 
would  lead  to  incorrect  results  in  CDI  calculations. 

To  explore  these  questions,  an  approximate  theoretical  model 
of  the  TI  wide  aperture  parabolic  localizer  system  was  developed. 
This  model  was  then  applied  to  the  prediction  of  CDI  for  a  poten¬ 
tial  installation  at  San  Francisco  International  Airport  and 
compared  to  predictions  for  the  same  installation  based  on  measured 
patterns . 

4.2  THE  TEXAS  INSTRUMENTS  PARABOLIC  ANTENNA  SYSTEM 

Figure  4-1  illustrates  the  positioning  of  the  components  of 
the  TI  localizer  system.  The  radiating  system  consists  of  seven 
directional  radiating  elements  and  a  parabolic  reflecting  screen. 
Three  of  the  elements  form  the  "directive"  array  and  radiate 


Figure  4-1.  Sketch  of  Parabolic  Localizer  Layout 

primarily  by  reflection  from  the  screen.  The  remaining  four 
elements  are  arranged  in  a  linear  array  and  are  oriented  so  as  to 
radiate  directly  to  the  landing  approach  path.  The  parabola  has 
a  focal  length  of  44.5  feet  (5A),  an  aperture  of  176  feet  (20A) 
and  a  height  of  18.5  feet  (~2A).  The  center  element  of  t.,e  direc¬ 
tive  array,  located  at  the  parabolic  focus,  generates  the  course 
carrier  plus  sidebands  signal.  The  course  sidebands-only  signal 
is  radiated  by  the  directive  array  side  elements,  which  are  placed 
slightly  in  back  of  the  focal  plane. 

A  detail  of  the  spacing  of  the  array  elements  is  shown  in 
Figure  4-2.  Each  element,  except  the  center  element  of  the 
directive  array,  consists  of  a  radiator  and  a  reflector.  This 
type  of  antenna  has  been  given  the  designation  "V- frame”  by  the 
manufacturer.  The  center  element  of  the  directive  array  is  of  the 
same  type  with  an  additional  component,  a  "director,"  to  give  it 
a  narrower  pattern.  The  pattern  of  the  V- frame  elements  of  the 
linear  array,  as  deduced  from  graphs  furnished  by  the  manufacturer. 
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Figure  4-2.  Detail  of  V- Frame  Element  Positioning 

is  tabulated  in  Table  4-1.  The  element  patterns  of  the  directive 
array  are  treated  differently  as  will  be  indicated  below. 

The  clearance  signal  is  generated  by  the  combination  of  the 
directive  and  linear  arrays.  The  contribution  of  the  center 
element  of  the  directive  array  may  be  adjusted  so  that  the  clear¬ 
ance  C+S  signal  has  a  near  null  on  the  course  line.  This  has  the 
virtue  of  enhancing  the  capture  effect  rejection  of  clearance 
multipath  interference. 

Patterns  claimed  by  the  manufacturer  for  a  system  so  adjusted 
are  shown  in  Figure  4-5.  To  achieve  these  patterns  the  signal 
amplitudes  fed  to  the  -..srious  elements  of  the  directive  and  linear 
arrays  are  adjusted  so  that: 

1.  The  on-course  strength  of  the  linear  array  C+S  radiation 
is  10  dB  below  the  course  C+3  signal  from  the  directive 
array. 

2.  The  clearance  C+S  of  the  directive  array  is  adjusted  to 
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TABLE  4*1  V- FRAME  HORIZONT. 


i  II ATI  ON  PATTERN 


Angle 

Relative 

Amplitude 

*  gle 

Relative 

Amplitude 

0 

1.000 

70 

.695 

±10 

.998 

80 

.  650 

±20 

.976 

90 

.588 

50 

.954 

100 

.  559 

40 

.912 

110 

.549 

50 

.842 

120 

.524 

60 

.772 

150 

.417 

- 

- 

180 

.435 

cancel  the  linear  array  C+S  signal  on  course. 

5.  The  relative  amplitudes  fed  the  four  elements  of  the 
linear  array  are: 

C+S:  0.0  .35/0°  .53/0°  0.0 

SO:  .165/270°  .55/270°  .55/90°  .165/90° 

4.  The  clearance  SO  signal  radiated  by  the  directive  array 
is  adjusted  to  equal  the  linear  array  SO  signal  at  5°. 

The  resultant  clearance  signal  radiation  patterns  are 
given  in  Figure  4-4. 

4.5  MATHEMATICAL  ANALYSIS  OF  REFLECTION  FROM  A  PARABOLIC  SCREEN 

The  radiation  contributions  of  the  four  element  linear  array 
can  be  calculated  m  the  usual  wav  by  standard  array  equations. 

In  order  to  generate  the  radiation  pattern  of  the  directive  array, 
we  require  approximate  formulas  for  diffraction  by  the  parabolic 
screen.  The  screen  is  formed  by  a  closely  spaced  horizontal  wire 
grid,  h'e  assume  that  it  may  be  treated  as  a  solid  perfect  reflector 


Figure  4-3.  Directive  and  Linear  Array  Clearance 
Radiation  Patterns 


Figure  4-4.  Combined  Clearance  Radiation  Pattern 
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to  the  horizontally  polarized  radiation  of  the  directive  array. 

The  elements  of  the  array  are  treated  as  point  source,  elevated 
9  feet  above  the  ground  plane  and  positioned  as  shown  in  F.gure 
4-2. 

To  develop  the  approximate  far-field  radiation  pattern  for  a 
point  source  reflected  by  the  parabola,  we  adopt  the  coordinate 
system  shown  in  Figure  4-5.  The  reference  origin  is  established  at 
the  point  where  the  parabola  focal  line  intercepts  the  ground  plane. 
From  this  point  the  vector  r'  gives  the  position  of  the  source 
element;  the  vector  r  gives  the  observation  point,  and  the  vector 
Tp  defines  a  variable  point  on  the  parabola  surface.  It  is  always 
considered  that  jr|  »  irp|  »  } r * | .  The  unit  normal  reference 
vector  on  the  parabola  n  is  taken  to  point  away  from  the  local 
center  of  curvature. 

The  reflected  electric  field  amplitude  E(r)  at  the  observation 

4 

point  is  taken  to  be  given  by  a  vector  Kirchhoff  integral: 

E(r)  =  ?  x  [f  x/fn  x  H.)  ds]  .  (4.1) 

Here  H^(rp)  is  the  magnetic  field  amplitude  incident  on  the  para¬ 
bola  due  to  the  source  element,  ip  is  the  Green's  Function: 


ikjr-rj  / 

*  =  e  /l?‘rpl  ’ 


(4.2) 


and  dS  is  the  increment  of  integration  area.  The  integral  in  (4.x) 
can  be  recast  using  the  radiation  field  relation: 


(4.5) 


where  n.  is  a  unit  vector  giving  the  direction  of  incidence  and 
E.(rp)  is  the  incident  electric  field  amplitude.  The  integral  then 
becomes 

E(r)  =  j  r  x  ^r  x  j  n  x(n^  x  E^)  -4>  dsj  .  (4.4) 

Ke  assume  that  the  incident  electric  field  is  horizontally  polarized 
and  may  be  expressed  in  the  general  form: 


where 


Ei(V  =  *s  A("i>  ”i  x  £  ’ 


iklr  -r'|  / 

*s  =  e  P  /iV 


(4 .  5) 


(4-6) 


A ( n . )  gives  the  field  amplitude  radiated  by  the  source  as  a  func¬ 
tion  of  direction.  Introduction  of  this  form  into  (4.4)  leads  to 
a  multiple-vector  product  which  may  be  reduced  under  the  assump¬ 
tions  that  n  is  horizontal  and  r  and  n^  are  nearly  so: 


(j  x(”  x(»i  x("i  x  *)))) 

A  ^  A  A  A  A  A  ^  A  A 

=  r  x  k  r*n  -  rxn^n^-kn-r 

A  A  /a  A  a  SH  A  a\ 

+  r  x  n  (nrk  nrr  -  r-kj 


=  r  x  k  r-n 


This  allows  (4.4)  to  be  written  as 


Efr)  =  |  r  x  k  J  d3  r-n  A(ni)  ps  # 


(4.7) 


(4-8) 


We  next  simplify  the  arguments  of  the  functions  $  and  i>,  i.e. 

|rp-r’{  and  jr-rpj,  recalling  the  assumption  r'  «  rp  «  r.  It  is 

convenient  to  put  rp  =  p  +  zpk,  with  p  the  horizontal  projection 

of  r  ,  and  then: 

P 


I  t_  -  r  *  I  -  p 


(-„  -  z'V 


!r  -  r  i  =r-r-p-z_r-k+- 


(4.9) 


(4.10) 


where  terms  such  as  jr’j~/p  have  been  rejected  as  negligible.  In 
order  to  factor  the  integrand  of  (4.8),  we  further  simplify  the 
product  $  $  after  inserting  the  approximations  above  and  we 

obtain: 
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(4-11) 


ikr 

e 


Negligible  terras  have  been  introduced  in  the  third  exponential 
factor  in  order  to  complete  the  square  appearing  there.  Since  the 
vector  p  describes  a  point  on  a  parabola,  it  satisfies  the  exact 


relations : 


(4.12) 


P 


+  2f 


(4.15) 


In  order  to  obtain  a  manageable  integral,  we  need  to  make  in 
certain  places  the  approximation  p  «  -  f.  The  element  of  integra¬ 
tion  combined  with  the  factor  r*n  can  be  reduced  in  the  following 
manner: 

A  A  ^  i  ^  a 

dS  r-n  =  k  •  dp  x  r  di 

p 

■  d»x  '  <f)x  d=,, 


=  (<f)v  ji  -  <f>x) 


do 


(4.14) 


Incorporating  these  results  into  the  integral  (4.S),  we  obtain  the 
form: 


E(f)  =  i 


ik(r  +  (l+ax)f) 


x  £  F, 


H 


rV 


(4 


15) 


where 


A(n.) 


ik 


d-ax) 


°v2 
47  ' 


v* 


(4.16) 
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dz?  e 


(4.17) 


a  =  r  ♦  r'/f 


=  (ax,  av,  a;) 


(4. IS) 


The  vertical  integral  F^,  is  a  Fresnel  integral  which  is  easily  put 
into  standard  fora  bv  introducing  the  variable 


u  =  \2/('a± }  (sp-s*-f  r-s)  - 


(4.19) 


However,  up  to  this  point,  this  integral  does  not  account  for  any 
reflections  fros  the  ground  plane.  Terms  aust  be  added  to  include 
the  image  of  the  source  in  the  incident  field  and  the  image  of  the 
screen  in  the  reflected  field.  With  these  terms  F^  can  be  reduced  to: 


r P  IS-U  -r'-f  J-i)2  1S(;  ‘i'-t  f-i)2 

PV  ■  i  4lp  *  '  C 

•"p 


(4-20) 


with 

Sy  =  k/2f  . 

The  horizontal  integral  fp  is  not  as  easily  treated-  We  note, 
however,  that  if  the  expression  1-a^  is  sufficiently  snail  the 
quadratic  ters  of  the  exponential  can  be  ignored-  This  appears  to 
be  the  case  when  the  deviation  angle  froa  the  course  to  the 
receiver  is  as  auch  as  14°.  This  is  s  large  enough  range  to  cover 
most  of  the  sector  in  which  the  radiation  of  the  directive  array 
is  significant.  Therefore,  we  proceed  by  dropping  the  quadratic 
t ere  of  the  exponent,  which  allows  Fp  to  be  treated  as  a  Fourier 
integral.  Putting  a  =  2f  and  assuming  an  expan; ion  of  A(n-): 


A(n.)  =  A  +  A, (o„/2£) ~  * 


(4-21) 
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Figure  4-6 I  Theoretical  Course  and  Clearance 
Signal  Patterns 


range  up  to  about  6°.  Such  differences  are  not  serious.  To 
confirm  this,  the  theoretically  calculated  parabolic  localizer 
patterns  were  inputted  to  the  ILSLOC  progran  as  measured  patterns 
with  no  assumed  phase  differences.  The  resulting  orbit  CDI  com¬ 
pared  with  the  model  generated  orbit  CDI  shows  no  significant  dif¬ 
ference  (Figure  4-7). 

All  of  this  suggests  that  the  TI  parabolic  localizer  signal 
patterns  can  be  adequately  generated  in  the  IL5LGC  program  from 
measured  data.  If  this  approach  is  taken  the  vertical  dependence 
of  the  simulated  radiation  will  correspond  to  a  discrete  antenna 
height.  Reduction  of  the  theoretical  vertical  dependence  (Equa¬ 
tion  (4.20))  in  the  limit  of  small  elevation  angles  indicates  that 
the  height  of  the  V- frame  elements  above  the  ground  plane  is  the 
correct  choice. 


5.  APPLICATIONS  OF  ILS  PERFORMANCE  PREDICTION  MODEL 


S.l  INTRODUCTION' 

The  Transportation  Systems  Center  ILS  local  iter  performance 
prediction  model  was  used  to  predict  ILS  localizer  performance  at 
several  different  airports.  Specifically,  the  sodei: 

1.  Nas  applied  to  Syracuse-Hancock  International  Airport  to 
obtain  sodei  validation 

2.  Nas  used  to  determine  which  type  of  localizer  was  needed 
to  produce  acceptable  course  structure  in  the  presence 
of  a  lar§e  747  hangar  located  near  runway  2S  at  the  San 
Francisco  Airport 

3.  Nas  used  to  predict  what  the  ILS  localizer  perforsance 
would  be  if  a  proposed  limited  access  roadway  were  to  be 
constructed  in  front  of  the  local Izer  antenna  at  the  New 
Orleans  Airport 

4.  has  used  to  determine  comparative  ILS  localizer  performance 
for  a  proposed  airport  layout  at  the  Bailas/Fcrt  North 
Airport  site. 

These  varied  applications  of  the  model  illustrate  the  poten¬ 
tial  usefulness  of  the  sodei: 

1.  As  an  aid  to  airport  planners  and  architects  by  being 
able  to  predict  the  effect  a  proposed  airport  layout  has 
on  ILS  performance  and  suggesting  possible  modifications 
(e.g.,  Dallas/Fort  North) 

2.  As  an  aid  to  particular  airport  managers  in  predicting 
the  effect  a  proposed  neu  hangar  or  roadway  has  on  ILS 
performance  and  suggesting  specific  location  or  orienta¬ 
tion  modifications  of  the  proposed  structure  for  accept¬ 
able  ILS  performance  (e.g..  New  Orleans) 

3.  As  an  aid  tc  the  FAA  regions  in  deciding  on  the  relative 
performance  of  different  ILS  localizer  antenna  candidates 
for  runwav  instrumentation  or  run wav  ungrading  (e.g.. 


San  Francisco  and  Syracuse  Airports) .  Each  of  these 

applications  of  the  model  will  now  be  discussed  (more 

extensive  information  may  be  found  in  individual  reports 

5-8 

previously  issued). 

5.2  INPUT  DATA  REQUIRED 

The  input  data  to  the  computer  program  for  these  applications 
follow  the  general  format  outlined  below,  (see  Reference  3  for  a 
complete  description  and  listing  of  the  computer  program). 

The  particular  localizer  antenna  desired  is  chosen  first.  ’.Ve 
have  the  choice  of  a  V-ring,  8-loep,  waveguide;  a  measured  pattern 
(single  or  capture-effect  system)  and  a  theoretical  pattern  (single 
or  capture  effect  system),  We  then  choose  the  transmitter  fre¬ 
quency  (usually  110  MHz).  We  place  the  chosen  antenna  at  the  ori¬ 
gin  and  input  the  distance  between  antenna  and  runway  threshold. 

This  is  used  for  both  flight  path  orientation  and  for  course  width 
determination.  The  antenna  height  above  the  ground  is  next  required 
as  well  as  the  flight  path  distance,  the  desired  spacing  between 
sample  points,  the  angle  of  approach  of  the  aircraft,  the  glide 
angle  and  the  height  of  the  aircraft  at  threshold.  For  an  orbit 
run,  the  starting  and  ending  angle  are  input  as  well  as  the  desired 
angular  spacing  between  sample  points  and  the  radius  and  height  of 
the  orbit.  The  velocity  of  the  aircraft  is  chosen  next  which  is 
used  to  calculate  the  Doppler  effect  on  the  receiver. 

After  inputting  the  antenna  patterns  used  at  the  location 
studied  as  well  as  the  assumed  trajectory  of  the  aircraft,  the 
location  and  orientation  of  the  derogating  jrfaces  are  next 
delineated.  The  shape  of  the  scatterer  is  input  first,  rectangular 
or  cylindrical,  followed  by  the  location  coordinates  of  the  scat¬ 
terer.  The  angle  made  by  the  base  of  the  object  relative  to  the 
runway  centerline  (orientation)  is  input  next,  followed  by  the  angle 
of  tilt  of  the  object  and  its  width  and  height.  Walls  that  are 
too  large  are  automatically  divided  into  smaller  ones  based  on  a 
minimum  error  criterion. 
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After  the  input  of  this  data,  the  electric  field  at  the  sur¬ 
face  of  the  scatterer  due  to  the  direct  and  ground  reflected  signals 
is  calculated.  Following  this,  the  electric  field  at  each  point 
along  the  flight  pzth,  due  to  the  direct,  ground  scattered  and 
building  scattered  signals,  is  calculated.  The  program  then  loops 
back  permitting  the  summation  of  the  effects  of  many  scatterers. 
Finally,  the  program  calculates  the  CDI  that  would  be  seen  at  each 
receiver  point. 

The  looping  permits  the  repetition  of  a  run  with  changes  in 
some  or  all  of  the  variables.  For  example,  a  run  may  be  repeated 
with  the  same  antenna  but  with  a  new  set  of  scatterers  or  with  the 
same  set  of  scatterers  but  with  a  different  antenna. 

A  second  (mini)  program  takes  the  output  tape  generated  by 
this  main  program  and  converts  it  to  the  observed  CDI  by  simulating 
the  effect  of  a  low  pass  filter  (to  so-called  "dynamic”  results  as 
opposed  to  the  "static"  or  instantaneous  CDI  values  obtained  as 
output  from  the  main  program).  The  receiver  time  constant  and  air¬ 
craft  speed  must  be  inputted  into  this  program. 

Finally,  a  plotting  routine  program  generates  a  variety  of 
graphs  of  the  static  and  dynamic  CDI's  as  well  as  of  the  antenna 
patterns. 

5.3  SYRACUSE -HANCOCK  INTERNATIONAL  AIRPORT 

In  the  case  of  Syracuse-Hancock  International  Airport,  flight 
test  recordings  were  available  against  which  the  model  could  be 
compared.  Syracuse  Airport  is  an  example  of  an  airport  with  a 
moderately  complex  environment  as  shown  in  Figure  5-1,  whose  in¬ 
strumented  runway  number  28  is  FAA  rated  for  Category  I  instru¬ 
mented  flights.  Pictured  in  Figure  5-1  are  a  number  of  structures 
near  runway  28,  numbered  from  1  to  18,  which  are  sources  of  scat¬ 
tered  ILS  localizer  radiation.  Since  the  localizer  radiation 
striking  and  being  reflected  from  these  structures  makes  an  angle 
other  than  zero  with  the  runway  ccnterlir.e,  the  radiation  striking 
an  aircraft  receiver  antenna  which  is  on  line  with  the  runway 
centerline  will  contain  unbalanced  150  Hz  and  90  Hz  modulations. 
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st rumen ted  Runway  28 


Thus,  the  aircraft  will  receive  indications  of  an  of f-Cv-nterline 
approach  due  to  the  combination  of  these  reflected  signals  with 
the  direct  on-course  signal.  In  order  to  approximate  the  amount 
of  derogation  caused  by  the  signals  reflected  from  the  buildings 
shown  in  Figure  5-1,  the  buildings  were  modeled  as  perfectly  re¬ 
flecting  rectangular  walls. 

The  localizer  is  an  eight-loop  antenna  located  2000  feet 
beyond  the  end  of  the  runway  at  an  effective  height  of  12  feet 
above  ground.  The  tailored  course  width  is  3.64°,  and  the  glide 
path  is  2.5°.  Kith  this  data,  the  couri>e  deviation  indication 
(GDI)  in  microamperes  was  calculated  and  compared  with  a  flight 
recording  taken  in  January,  1972.  This  is  shown  in  Figure  5-2. 

The  agreement  is  good,  differences  being  due  to  any  of  a  number 
of  causes: 

o  The  approximation  of  perfectly  reflecting  walls; 

o  The  assumption  that  the  aircraft  was  precisely  on  a  2-S° 
glide  path  and  was  travelling  at  120  knots  on  a  runway 
centerline  approach; 

o  The  assumption  that  all  building  doors  were  closed. 

The  CD I  flight  test  recordings  and  theoretical  predictions 
show  that  runway  28  meets  Category  I  FAA  requirements  but  not 
Category  II.  The  eight- loop  localizer  antenna  used  at  Syracuse 
Airport  has  a  rather  broad  azimuthal  beam  which  illuminates  many 
of  the  buildings  shown  in  Figure  5.1  causing  the  reflections 
responsible  for  the  indicated  course  deviation.  If  an  antenna 
whose  radiation  pattern  fell  off  more  rapidly  with  azimuthal 
angle  were  used,  several  of  the  buildings  which  arc  now  illumin¬ 
ated  by  the  eight-loop  might  fall  outside  the  range  of  the  nar¬ 
rower  azimuthal  beam  antenna  leading  to  a  better  course  structure 
for  the  runway. 

It  is  known  that  the  V-ring  antenna  pattern  energy  falls  off 
more  rapidly  with  azimuthal  angle  than  does  the  eight- loop  antenna 
The  V-ring  antenna  was  therefore  substituted  for  the  eight-loop 
in  the  theoretical  model  and  the  CDI  recalculated.  As  expected. 
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figure  5-2.  Comparison  Between  Flight  Recording  and  Theoretical  Prediction 
of  Course  Structure,  Runway  28,  Syracuse-Hancock  Airport 


the  recalculated  CDI  is  considerably  reduced  from  that  which  was 
obtained  with  the  eight- loop  as  shown  in  Figure  5-3.  In  fact, 
runway  28  at  Svracuse-Hancock  Airport  may  be  upgraded  from  an 
instrumented  Category  I  runway  to  one  which  meets  FAA  requirements 
for  Category  II  by  replacing  the  existing  eight-loop  localizer  with 
a  V-ring  localizer. 

5.4  DALLAS /FORT  WORTH  STUDY 

In  the  case  of  Dallas/Forth  Worth  Airport,  input  data  for  the 
model  were  obtained  from  a  set  of  architectural  drawings  of  the 
proposed  four  instrumented  runways  and  Phase  1  building  construc¬ 
tion  plans.  This  data  was  used  to  predict  and  compare  ILS  local¬ 
izer  performance  using  capture  effect  and  non-capture  effect 
localizer  antenna  systems  for  the  four  runways.  The  proposed  air¬ 
port  layout  plan  is  shown  in  Figure  5-4.  Runways  17R  and  35L 
running  north  and  south  are  located  to  the  west  of  the  north-south 
(eighteen  lane)  highway  shown  flanked  by  terminal  buildings  (semi¬ 
circles  in  Figure  5-4)  located  in  the  areas  marked  as  OW,  OH,  JW,  IE, 
etc.  In  the  Phase  1  building  plan,  only  terminal  buildings  in  the 
areas  marked  2E,  2W,  3E  and  4E  will  be  constructed.  In  addition 
to  the  construction  of  these  large  3000-foot  terminal  buildings,  a 
124-foot  high,  320-  by  60-foot  hotel  will  be  constructed  in  loca¬ 
tion  3W,  as  well  as  a  175 -foot  high  tower  located  symmetrically 
between  the  17R/35L  and  17L/55R  runways  and  a  Braniff  and  a  Delta 
hanger  located  as  shown  in  the  figure.  Runways  17I./55R  also  run 
north  and  south  and  are  located  to  the  east  of  the  highway.  The 
x-raarks  in  the  figure  represent  the  proposed  locations  of  the 
localizer  antennas.  The  x-raark  circumscribed  by  a  circle  repre¬ 
sents  a  proposed  Category  II  installation  serving  runway  17L. 

The  performance  of  a  V-ring  localizer  and  an  Alford  IB  local¬ 
izer  were  compared  for  the  four  runways.  The  CDI  was  predicted 
along  a  centerline  aircraft  approach  at  a  2.5°  glide  path.  A 
representative  result  comparing  an  Alford  IB  and  a  V-ring  localizer 
for  runway  35L  is  shown  in  Figures  5-5  and  5-6  (the  complete  study 
may  be  found  in  Reference  5).  These  two  figures  show  the  dynamic 
CDI’s  (assuming  a  receiver  time  constant  of  0.4  and  an  aircraft 
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Theoretical  CBI's  at  Syracuse-Hancock  Airport 
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Theoretical  CDI*s  at  Syracuse-Hancock  Airport 
Using  a  V-Ring  Localizer 


Figure  5-5 


Coizparison  Between  Course  Structures  Using 
a  V-Ring  and  8-Loop  Localizer  Antenna 


Proposed  Dnllas/Port  Worth  Airport  Layout  Plan 


Figure  5-5.  Predicted  Course  Structure  in  Presence  of  Terminals, 
Hotel  and  Tower  Using  a  V-Ring  Array  for  Runway  SSL, 
Dallas/Port  Worth  Airport 
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Plguro  S-6.  Predicted  Course  Structure  in  Presence  of  Terminals,  Hotel 
und  Tower  Using  an  Alford  111  Loco  User 


speed  of  120  knots)  predicted  by  the  model  in  the  presence  of 
terminals  2E,  2K,  5E  and  4E,  the  hotel  and  the  tower.  Additional 
derogation  due  to  the  Braniff  hangar  when  using  the  Alford  IB  array 
is  shown  in  Figure  5-7.  The  derogation  due  to  the  hangar  is  snail. 
(The  information  on  the  Delta  and  Braniff  hangars  case  too  late  to 
include  in  the  original  CDI  prediction  containing  the  terminal 
buildings,  hotel  and  tower  shown  in  Figures  5-5  and  5-6.  At  worst, 
the  effect  of  the  two  hangars  (only  the  Braniff  hangar  is  shown 
here)  will  be  to  add  to  the  total  CDI  shown  in  Figures  5-5  and 
5-6,  and  using  an  Alford  IB,  their  effect  should  be  negligible.) 

The  next  two  figures,  5-8  and  5-9,  show  the  orbit  clearance 
for  runway  3SL  taken  for  an  aircraft  flying  600  feet  above  the 
ground  and  25,00G  feet  from  the  localizer.  Both  the  V'-Ring  and 
Alford  array  patterns  meet  the  FAA  specifications,  yielding  a  CDI 
of  150  microamperes  or  sore  within  the  azimuthal  sectors  +35°. 

Ke  concluded  from  the  application  of  the  model  to  the  Dallas/ 
Fort  iforth  Airport  plan  that  the  capture  effect  system  tAlford  IB) 
should  be  used  in  preference  to  the  V-Ring  to  insure  that  FAA 
Category  I  and  Category  II  requirements  be  met  on  their  respective 
runways  under  ifce  Phase  1  construction  plan. 

5.5  SAN  FRANCISCO  AIRPORT  STUDIES 

The  Transportation  Systems  Center  model  was  also  applied  to 
runway  2S  at  the  San  Francisco  Airport  in  order  to  help  determine 
an  acceptable  localizer  antenna  tor  the  course-  Large  derogation 
to  the  ILS  localizer  performance  was  expected  here  because  of  the 
presence  of  a  large  747  American  Airlines  hangar  located  relatively 
close  to  the  runway.  Several  different  localizer  systems  were 
tried:  A  V-ring,  the  Alford  IB  or  14/6  array,  an  Alford  2A  or 
22/S  array  and  the  Texas  Instrument  Parabolic  localizer  antenna. 

In  addition,  the  Alford  and  Texas  Instruments  antennas  were  oper¬ 
ated  at  different  capture  effect  ratios  than  is  normally  done  in 
order  to  improve  their  performance. 

The  scattering  situation  is  depicted  in  Figure  5-10.  The 
hangar  size  is -560  feet  by  456  feet,  135  feet  high  and  is  oriented 
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Figure  5-9-  Predicted  Orbit  Clearance  Structure  in 
Presence  of  Terainals,  Hotel  and  To wer 
Using  an  Alford  IB  Localizer 
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Figure  5-10.  Schematic  of  Scattering  Situation  for 
Runway  2$,  San  Francisco  Airport 


parallel  to  and  1328  feet  away  from  the  runway  centerline.  The 
CD1  was  predicted  for  an  aircraft  flying  a  2.5°  glide  slope  center- 
line  approach. 

Figure  5-11  shows  the  predicted  CDI  when  a  V-ring  is  used. 

The  course  is  very  poor  with  this  antenna  owing  to  the  large  amount 

of  sideband  energy  reflected  from  the  hangar  into  the  flight  path 

of  the  aircraft.  The  reason  for  this  is  because  of  the  relatively 
large  V-ring  sideband  radiation  energy  existing  at  13°,  the  approx¬ 
imate  angular  position  of  the  hangar  off  centerline. 

Figures  5-12,  5-13  and  5-14  show  the  predicted  CDI  when  using, 

in  the  normal  mode  of  operation,  an  Alford  14/6,  an  Alford  22/8 

and  the  Texas  Instruments  localizer  antenna,  respectively.  The 
course  structure  has  improved  significantly  in  using  these  capture 
effect  systems,  but  still  remains  poor.  In  order  to  further  im¬ 
prove  the  course  structure,  we  operated  the  Alford  22/8  array  and 
the  Texas  Instruments  Parabolic  localizer  at  a  higher  course  to 
clearance  capture  effect  ratios  than  is  normally  done.  We  raised 
the  capture  effect  rati  from  10  db  of  normal  operation  to  16  db. 

The  results  are  shown  4  •»  Figures  5-15  and  5-16.  The  results  show 
that  significant  improvement  is  attained  in  this  higher  capture 
effect  ratio  mode  of  operation.  A  localizer  performance  is  pre¬ 
dicted  which  allows  FAA  Category  I  requirements  to  be  met.  and, 
marginally,  also  Category  II  requirements.  However,  we  caution 
that  it  has  not  been  established  that  either  the  Texas  Instruments 
Parabolic  antenna  or  the  Alford  22/8  array  can  function  reliably 
at  the  higher  16  db  capture  effect  ratio  used  in  our  theoretical 
model.  Whether  or  not  such  operation  is  possible  has  to  be 
determined  in  the  light  of  the  power  output  limitations  of  the 
course  and  clearance  transmitters. 

5.6  NEW  ORLEANS  AIRPORT  STUDY 

In  the  case  of  New  Orleans  Airport,  we  were  asked  to  predict 
the  performance  of  an  Alford  IB  localizer  for  runway  10  under  the 
assumption  that  a  roadway  would  be  built  in  front  of  this  localizer. 
The  roadway  would  be  a  limited  access  road  containing  automobile 
traffic  which  could  conceivably  derogate  localizer  performance. 
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V-Ring:  Dynamic  Flyability  Run 


Course  Structure  Using 
San  Francisco  Airport 
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Figure  5-14.  Predicted  Course  Structure  Using  the 
Texas  Instrument’s  Parabolic  Antenna 
Operating  with  a  10  dB  Capture  Effect 
Ratio,  Runway  28,  San  Francisco  Airport 
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Figure  5-15.  Predicted  Course  Structure  Using  the  Alford  2A  Localizer, 
Operating  with  a  16  dB  Capture  l- f feet  Ratio 
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Figure  5-16. 


Predicted  Course  Structure  Using  the 
Texas  Instrucents  Parabolic  Antenna 
Operating  with  a  16  dB  Capture  Effect 
Patio 
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Thr^i  =?*«rnate  installations  were  proposed  for  the  roadway,  one 
in  which  the  road  was  level  with  the  ground,  and  two  in  which 
the  road  would  be  located  at  different  depths  below  ground  level. 
Only  the  level  roadway  was  considered  in  this  Study- 


Figure  5-17  shows  a  partial  airport  layout  with  the  proposed 
United  access  road  running  in  a  north-south  direction  perpendicu¬ 
lar  to  the  instrumented  runway  number  10.  The  local  iter  sketched 
in  the  figure  is  to  be  an  Alford  IB  array  located  850  feet  from 
the  roadway  centerline  at  its  closest  point.  Various  types  oi 
automobile  traffic  and  one  case  of  a  single  trailer  truck  were 
assumed  and  the  resulting  Alford  IB  local  iter  operating  with  a 
3.59'*  course  width  and  a  10  dB  capture  effect  ratio,  was  predicted 
for  an  aircraft  flying  a  centerline,  2.5°  glide  slope  approach. 

The  types  of  vehicular  traffic  for  which  the  local  iter  per¬ 
formance  was  predicted  consisted  of  a  single  perfectly-reflecting 
two-sided  car,  20  feet  by  8  feet  and  6  feet  high  located  on  the 

O 

road  so  as  to  give  a  worse  case  single- car  result. 

A  single  perfectly-reflecting  two-sided  trailer  truck  45  feet 
by  S  feet  and  15  feet  high  was  also  modeled  for  comparison  pur¬ 
poses  with  the  car.S 

A  traffic  jam  along  part  of  the  roadway  was  also  considered. 

The  traffic  jam  situation  was  modeled  as  a  long  perfectly- reflect- 

8 

ing  wall  6  feet  high. 

Finally,  a  worse  case  car  spacing  (called  "Fresnel  cars") 
was  considered.  Each  car  in  this  traffic  spatial  distribution  is 
assumed  to  have  two  perfectly-reflecting  sides  20  feet  by  S  feet 
and  6  feet  high,  located  as  shown  in  Figure  5-lSa. 

The  CDI's  for  these  vehicular  traffic  situations  were  calcu¬ 
lated  along  the  flight  path  from  a  distance  of  25,000  feet  from 
the  localizer  to  within  2000  feet  of  it.  In  all  cases  the  deroga¬ 
tion  due  to  the  presence  of  the  representative  traffic  distribu¬ 
tions  modeled,  is  very  small,  indicating  that  the  proposed  (level) 
roadway  carrying  predominately  automobile  traffic  will  add  only 
insignificant  amounts  of  additional  derogation  to  the  course 
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structure  for  runway  10  using  an  Alford  IB  localizer.  The  predicted 
GDI  in  the  presence  of  the  "Fresnel  Car  Spacing"  is  shown  in  Figure 
5-lSb. 

g 

As  noted  in  the  report  on  the  Sew  Orleans  study,  care  must  be 
taken  in  accepting  these  results  without  additional  conf irs3t ion. 
This  is  because  the  TSC  model  was  developed  for  the  prediction  of 
ILS  perforaance  when  the  main  derogation  is  caused  by  scattering 
fros  large  objects  suen  as  hangars  where  diffraction  effects  are 
of  secondary  importance.  In  the  present  problem,  ILS  derogation 
is  due  to  scattering  fros  small  objects  (snail  cospared  to  several 
wavelengths)  where  diffraction  effects  are  important-  It  is  there¬ 
fore  difficult  to  predict  how  accurate  the  TSC  aodel  is  in  this 
case  and  independent  ‘validation  is  necessary. 

5-7  SUMMARY  OF  APPLICATIONS 

In  summary,  the  examples  cf  the  applications  of  the  model 
cited  indicate  the  potential  usefulness  cf  the  aodel. 

•  in  a  relatively  inexpensive  way  we  were  cble  to  suggest 
that  a  Category  I«  operation  of  the  instrumented  runway 
at  Svracuse-Hancock  Airport  Is  possible  with  the  replace- 
sent  of  the  present  eight- loop  localizer  antenna  array  by  a 
V-ring  antenna  array. 

•  for  the  Dallas/Fort  Worth  Airport  we  showed  that  the  pro¬ 
posed  Phase  1  construction  plan  was  commensurate  with 
acceptable  ILS  localizer  performance  when  a  capture  effect 
systea  (Alford  IB)  localizer  antenna  array  was  used. 

•  For  runway  28  at  the  S3n  Francisco  Airport,  we  made  a 
comparative  study  of  different  localizer  antenna  systems 
and, in  a  conservative  estimate  assuming  that  the  derogator 
was  perfectly  reflecting,  we  showed  that  Category  I  FAA 
requirements  would  be  met  only  by  a  capture  effect  system 
operating  at  a  higher  capture  effect  ratio  than  is  norm¬ 
ally  done. 
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•  Finally,  for  New  Orleans  Airport  we  suggested  that  the 

proposed  roadway  could  be  built  without  causing  significant 
derogation  to  the  operation  of  a  capture  effect  localizer 
system  though  additional  confirmation  was  necessary  in 
this  case. 


APPENDIX  A 


SCATTERING  FROM  A  PERFECTLY  CONDUCTING  VERTICAL 
CIRCULAR  CYLINDRICAL  SECTION  ELEVATED  ABOVE  GROUND 
In  Section  2.3,  localizer  signal  scattering  by  a  right,  circu¬ 
lar  cylinder  situated  on  the  ground  was  investigated.  In  this 
Appendix  the  scattering  formulas  appropriate  to  a  cylinder  which 
is  elevated  above  ground  level  are  derived.  Figure  A-l  shows  the 
geosetrv  for  this  scattering  calculation.  The  reference  point,  0, 
is  on  the  midpoint  of  the  axis  of  the  cylinder  with  section  length 
&z.  z^  is  the  height  of  the  reference  point  above  a  flat  infinitely 
conducting  ground  p’ane.  P  is  a  point  on  the  surface  of  the  cylin¬ 
der  and  r  is  the  vector  fros  0  to  P.  It  is  convenient  to  set  up 
the  coordinate  system  x*,  y',  z’  on  the  cylinder  such  that  y*  is 
an  extension  of  (the  ground  projection  of  R1 )  and  x*  is  per¬ 
pendicular  to  (see  Figure  A-2a).  As  shown  in  Reference  1  for 

saall  elevation  angles  of  the  rays  to  the  scattering  surface,  the 
incident  electric  field  is  given  by: 


Figure  A-l.  Geoaetry  for  Scattering  from 
Vertical  Circular  Cylinder 
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Hi  =  U  Eq  £(*’) 


2ikz *H/D  . 
1  -  e  pi 


(A-  1 J 


Where  U  is  a  unit  vector  giving  the  direction  of  electric  field 
polarization  and  f(i|>')  is  the  horizontal  antenna  gain.  The  angle 

•f 

<P'  is  the  angle  between  the  ground  projection  of  R  and  the  runway 
'•p '  is  approximated  by  ip  the  angle  between  the  ground  projection 
£  and  the  runway).  The  variable  z'  is  the  distance  of  the 
field  point  P  above  the  ground  plane  and  may  be  expressed  as: 

z'  =  zi  +  zc 

where  zc  is  the  distance  of  the  field  point  P  above  a  horizontal 
plane  cutting  the  center  of  the  cylinder  as  shown  in  Figure  A-2b. 
The  quantity  kR  may  be  expanded  as: 


kR  =  kRj  +  kR1 


.  t[r2  - 

•  x*  +  — — - - - A  <’  ♦ 


(A-  2) 


where  k  is  2rr  divided  by  the  wavelength  A.  Then,  under  the 
restriction  that: 


?  1 
kr  _  irr" 


21^  ~  isy 


(A-3) 


ive  will  approximate  kR  by  the  first  two  terms  on  the  right  hand 
side  of  Equation  (A-2).  In  accordance  with  the  geometry  shown  in 
Figures  A-l  and  A-2  the  following  equation  is  obtained: 


a  sin  4> *  cos  0j 


sin  0. 
c  1 


(A- 4) 


where  is  the  angle  between  x'  and  the  projection  of  r  onto  the 

horizontal  plane  through  the  reference  point,  and  where  0.  is  the 
-►  ■+  ■* 
angle  between  R^  and  the  projection  of  R^  on  the  ground  plane. 

The  electric  field  incident  at  point  P  may  then  be  written 
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The  incident  magnetic  field  is  given  by: 

H-  =  |E.|  /eTTT  (R  x  0) 


(A-  6) 


where  /p/e  is  the  wave  impedance  of  the  medium.  We  assume  U  to  be 
a  horizontal  vector  and  thus  obtain: 


H.  = 
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where: 


sin  9j  =  (Zj-HJ/Rj 


cos  9X  =  Dpl/R1 
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(A- 7b) 

(A-S) 
(A- 9) 

(A- 10) 


and  I)  is  the  ground  projection  of  R  and  e  ,  e  and  e_  are  the  unit 
P  x  y  £. 

vectors  in  the  x,  y,  and  z  directions,  respectively.  In  all  sub¬ 
sequent  ca’culations  the  small  elevation  angle  approximation  is 
made  so  that  the  second  terra  on  the  right  hand  side  of  Equation 
(A- 7b)  may  be  neglected. 


The  surface  current  is  given  bv: 


J  “  "total  x  n 


(A- 11) 


where  n  is  the  unit  normal  vector  pointing  into  the  surface.  To 
calculate  the  scattered  f  «*  field  we  require  the  Schelkunoff 
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radiation  vector: 
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(A- 12) 


where  the  integration  is  over  the  cylinder  surface.  When  the 
diameter  is  much  larger  then  the  wavelength  of  the  radiation  we  can 
use  Kirchoff’s  approximation  in  which  the  tangential  component  of 
the  total  magnetic  field  is  approximated  by  twice  the  tangential 
component  of  the  incident  magnetic  field  on  the  illuminated  side 
and  the  very  small  field  on  the  shadow  side  is  approximated  by  a 
null  field. 

Then  for  the  integration  on  the  illuminated  side  we  require 
the  following: 


J  5  2H  •  e.  x  n 

Z 1  i. 


(A- 15) 


where  H  .  is  the  vertical  component  of  the  incident  magnetic  field 
as  given  by  Equation  (A-7).  The  x1  and  y*  components  of  J  are 
then  (see  Figures  A- 2a  and  A- 2b) 


J  S  -  2H  •  sin 

X  ~  I 


J  =  -  2H  •  cos  <?’ 
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(A- 14a) 
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where  we  have  dropped  the  primes  on  the  x  and  y  subscripts  in  the 
expressions  (A- 14a)  and  (A- 14b).  We  also  have  the  relation: 


R7  =  l  cos  e7  sin  S  -  n  cos  g7  cos  B  ♦  e„  sin  e7 

where  X,  n  and  e_  are  the  unit  vectors  in  the  x*t  y*  and 
tions,  respectively,  and  where: 

sir  c,  = 
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Since: 


a  cos  £  -  a  sin  <J>*  n 


(A- 17) 


we  have 


R7  *  r  =  a  cos  c7  sin  8  cos  ♦  a  cos  c7  cos  S  sin 


+  2c  sin  e7 
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The  components  of  the  Schelkunoff  radiation  vector  in  x*  and 
y*  directions  are  given  by: 
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From  Equations  (A-5),  (A-7),  (A-1S),  and  (A-19)  the  following  are 
calculated: 


Y  =  Nx  cos  8  +  N  sin  8 
=  ay  I  IA 

O  y 


(A- 20a) 
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where 


r—  ikR,  / 

‘  2  Vj  61  Eo  f<*>  e  /»!  • 


(A* 21) 


SO 


I  =  Az  |sinc^k(sin8 ^-sinc?)Az/2j 


UkHZj/D 


pl  sine  [k(sin01-sine?+2H/Dpl)Az/2j |  ,  (A-22) 


where 


sine  x  =  sin  x/x. 


(A- 25) 


J 

i $  =  f  d$’  sin(0*+S)  e 


i(a  sin  O'  +  8  cos  O') 


,  (A- 24a) 


with 


a  -  -  ka(cos  6j  +  cos  t 2  cos  S) 


(A- 24b) 


SQ  =  -  ka  cos  e2  sin  6  . 


(A-24C) 


Using  the  identities  to  define  C  and  G: 


aQ  sin  O’  +  S0  cos  O'  =  C  cos  (O’-G)  , 


we  may  express  1^  as 


I,  =  cos  B  L  ♦  sin  S  I~ 

V  O  / 


(A -2 5a) 


=  C  sin  (A- 25b) 


(A* 26) 


Dropping  the  prime  on  we  have 


*  r 

I&  =  do  sin  0  cos^C  sin^+G-o^j  ♦  i  do  sin  0  sin^C  cosfc-Gij 


(A- 27) 


Using  the  expansion  (Reference  9) 
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(A- 28) 


cos  (C  sin  6)  =  JQ(C)  +  2  ^  J,K  (C)  cos  (2K6), 

K=1 

where  J  is  the  Bessel  function  of  1st  kind  of  order  n,  we  obtain 
n 

for  the  real  part  of  I,: 


Ren6)  J0(c)  /  d$  sin  <?  + 


I 


is 

2  J  dd  sin  e  cos(2K4>)J?K(C}  . 


sN  * c)] 

2  j  d$  sin  $  sin(2K<>) J?K(C)  sin|^2K^j-  ♦  g| 


from  which  is  obtained: 

CD 

W  '  2  VC>  *  I2 

K=1 


1  *  (-1) 


1+2K 


1  - 


- T 

(  2K)  “ 


J?K(C)  cos 


- 

2K0  *  Gj  (A- 29) 


Likewise  using  the  expansion: 


JinfC  cos  S)  =  2  ]>  (-)K  J?K+lfC)  cos  [(2K  +  U®]  (A-30J 


K=0 


we  obtain  for  the  imaginaiy  part  of  1^: 


IM(I6)  =  Jj(C)  n  sin  G 


For  I7  in  Equation  (A- 26)  we  obtain 


I  -  -  J  d£  cos  p  cos 

J0 


(A- 31) 


[C  s  in^*-G-^j|  +  i  di  cos  >  sin|c  cos(#  -  G)1 


(A- 32) 


Again  using  Equations  (A-28)  and  (A-50)  we  obtain: 


**  H 

Re(I?)  =  JQ(C)  j  d$  cos  %  *  2j7J.(C)  f  d<J  cos  $  cos(2K$). 
”0  0 

‘[“(I  *  GJ 


CCS 


2  J7K<C>  f  d*  cos  6  sinf2Kv)  sin[2K(f  +  G] 


ihe  first  two  terras  are  zero  so  that: 


R  (I,) 


- 2 1  j2k<c>-{  *  G!  (A-33) 


,  2K- 1* 


K=1 


(2K) 


and  similar  calculations  yield: 


Inl^  =  Jj(C)  -  cos  G 


(A-54) 


According  to  Equation  (A-2Sa)  where  the  constants  C  and  G  were 
introduced,  we  have: 


7  7\l/2 


C  '(“o  *  3o) 


(A-55a) 


=  ka 


cos2  Sj 


♦  2  cos  5,  cos  e7  cos  S  ♦  cos"  c7 


1/2 


sin  G  =  '^0/C 


cos  G  —  3  /C 
o 


(A- 55b) 
(A- 56a) 

(A- 56b) 


where  3  and  a  are  given  by  Equations  (A-24b)  and  (A- 24c),  and 
o  o 


G  =  Ira  G  ♦  i  sin  G)j 


(A -57) 
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where  R^  is  the  distance  from  the  reference  point  on  the  cylinder 
to  the  image  of  the  receiver  point. 


4  '  [Dp2  *  <z2*sl>2 


and  where 


(A-44) 


sin  el  =  (c-,+2  )/R; 


(A-45) 


And  I.  =  I  (see  Equation  A-22)  when  sin  el  replaces  -  sin  e,  there. 


Also  I.  =  I-  in  Equations  (A-24)  when 

if  0 


cos  4  =  Dp2/R2 


(A-46) 


replaces  cos  c?  there.  Finally,  yc  is  given  in  Equation  (A-21). 

The  total  induced  voltage  is  given  by  the  sue  of  the  contri¬ 
butions  froa  the  cylinder  and  its  isage: 


V  =  V  ♦  V 

'total  via 


(A-47) 


The  scattering  froa  the  complete  elevated  cylinderical  object 
is  then  obtained  by  siaply  sunning  up  the  contributions  from  each 
of  the  cylindrical  sections  treated  in  this  appendix. 
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APPENDIX  B 


(1)  SCATTERING  FROM  A  SLANTED  PERFECTLY 
CONDUCTING  TRIANGULAR  SLAB  ELEVATED 
ABOVE  GROUND 

(2)  SCATTERING  FROM  A  SLANTED  PERFECTLY 
CONDUCTING  DISK  ELEVATED  ABOVE 
GROUND 

B.l  SCATTERING  FROM  A  SLANTED  TRIANGLE 

In  a  previous  report.  Reference  1,  Section  4-4,  a  general 
formulation  for  the  scattering  from  a  perfectly-conducting  rectan¬ 
gular  flat  surface  was  developed.  That  formulation  is  modified  to 
obtain  the  scattering  from  an  elevated,  slanted  right  triangular 
slab  having  one  edge  parallel  to  the  ground  plane-  The  slab  is 
placed  above  a  flat  perfectly-conducting  ground  plane.  In  Equation 
(4.9?)  of  Reference  1  the  voltage  induced  in  a  horizontal  loop 
3ntenna  at  the  observation  point  by  the  scattered  field  from  a 
flat  slab  is  given  by: 


c  I,  cos 
o  b 


c7|cos  sin  c  cos  g  -  sin  cos  e  cos(S+y)J  (B-l) 


(Note  that  the  minus  sign  on  the  right  hand  side  of  this  equation 
was  inadvertently  omitted  in  Equation  (4.97)  of  Reference  1). 
Equations  (4.S2),  (4.83)  and  (4.86)  of  Reference  1  give: 


ikEof(f) 


exp  Jik(R2-*-R1)j  Ic/R2?-3 


(B-2) 


Ic  *  |exp|^ik(R1-r)J  ^1  -  exp(ik  2z’  H/D?1] 

•  exp(-ikRj-r)| 


ic  =  j&n  exp[ik(A-no)n]  exp[ik(8-£ou] 

•  exp  T  Ji!.  exp  ik(s  -  5o  -  -ip^S-£)  «]} 


exp  t  /  exp 


(B-5) 


(B-3a) 


(B-5b) 


. li*iili'»Wiii''|i1'lliiU'|l*|i|im 


(the  parameters  in  Equation  (B-5b)  are  given  in  Equations  (B-5) 
through  (B- 13)  below) . 


For  the  case  of  the  slanted  triangle  we  specify  the  boundaries 
of  the  surface  integration  in  Equation  (B-3).  Figure  B-l  shows 
the  general  geometric  relations  for  the  problem.  Figure  B-2a, 
b,  c  and  d  show  four  orien'ations  for  the  slanted  right  triangle. 
The  reference  point,  0,  is  taken  to  be  on  the  midpoint  of  the  side 
that  is  parallel  to  the  ground  plane.  For  all  four  cases  the  r, 
integration  goes  from  -  y  to  + 


L/2 

dn 

L/2  . 


CASE  1.  Acute  angles  pointing  to  the  left  and  upwards  (as 
shown  in  Figure  B-2a).  In  the  right  hand  side  of  Equation  (B-5) 
the  parameter  £  is  integrated  from  0  to  h(n/L  *  1/2).  W'e 
obtain: 


f  ikB"h/2 

Ic  ■  ( i/k)  <-  (1/B*)  e  L  sine  (K  L/2) 

+  (1/B")  cT  elkB"h/2  L  sine  (K’L/2) 

+  L  sine  [k(A-no)L/2][|r  -  eT]|  ,  (B-4) 

where: 


B"  =  B  -  £0  ♦  2H  sine  /DpJ  (B-5) 
B*  =  B  -  (B-6) 
B  =  cos  6^  cos  c  cos  t  +  sin  3^  sin  c  (B-7) 
£  =  -  cos  c j  cos  5  cos  £  ♦  sin  c,  sin  c  (B-S) 
nQ  =  cos  c7  sin  S  (B-9) 
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a.  Slanted  Right  Triangle 

b-  Slanted  Right  Triangl 

Pointing  to  the  Left 

Pointing  to  the  Right 

and  Upward 

and  Downward 

f 

4 

1 

i  juissi  ressr 

■  _ 

c.  Slanted  Right  Triangie 
Pointing  to  the  Right 
and  Upward 


d.  Slanted  Right  Triangle 
Pointing  to  the  Left 
and  Do wnwa rd 


Figure  B-2.  Four  Orientations  for  the  Slanted  Right 
Triangle 


CASE  4.  Acute  angles  pointing  to  the  left  and  downward  (Fig¬ 
ure  B-?,d).  The  parameter  £  is  integrated  from  -  h(n/L  -  1/2)  to  0. 
We  obtain: 


Ic  =  (i/k) 


{(i/B* 


)  e‘ik8M,/2  i  Sinc  (Y  L/2) 


-  (1/B")  eT'ikB"h/2  l  sinc  (Y’L/2) 


+  L  sinc 


[k(A-V>L/2](^  *  h  eT)}' 


(B-19) 


As  in  the  case  of  the  slanted  rectangular  scattering  formula¬ 
tion,  the  present  formulation  for  the  triangle  is  also  restricted 
to  the  case  when: 


riq  and  H“  <<:  1 


(B-20) 


where  X  is  the  radiation  wavelength. 

B.2  SCATTERING  FROM  AN  ELEVATED  SLANTED  PERFECTLY  CONDUCTING  DISK 

In  this  problem  we  take  the  reference  point  0  (see  Figure  B-3) 
to  be  at  the  center  of  the  disk.  The  boundary  of  the  disk  £'(n)  of 
radius  a,  is  given  by: 


2  2  2 

5’ (nr  +  n  = 


(B-21a) 


S*(n)  =  ±  Va2  -  n2 


(B- 21b) 


From  Equation  (B-3b)  we  have: 


a 

i*’j 


dn  exp(ikA')  - 


va  -  n 

/ 

J.n  ? 


12  2 

/a  -  n 


d£  exp(ikB’S) 
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exp  T 


/ 


dC  exp(ikB”C) 


■tr:  n2 


where  A'  =  A  -  n 


(B-22a) 


(B-22b) 


Integrating  over  C  we  obtain: 


where 


and 


lc  =  W 


21  ,  2  expT  I 


kB 


■/ 


I  a  ~  I  dn  f (n)  , 


(B-25) 


(B-24a) 


with 


ie  »  [  dn  F(n) 
-a 


f(n)  =  exp(ikA’n)  sin 


(kB'Va2  -  n2)  , 


(B-24b) 


(B-25a) 

F(n)  =  exp(ikA'n)  sin^kB"  Va2  -  n2^  (B-25b) 

The  integrals  1^  and  I£  may  be  approximated  by  the  following  scheme 


8 

5 

i=l 


Id  -  1  fC'.i)a/4 


(B-26a) 


>c  -  I  j)a/4 


(B-26b) 


i=l 
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Figure  B-3.  Diagram  Showing  Coordinate  Origin  of  £  and  n  Axis  a 
Center  of  Slanted  Disk.  Orientation  n  and  |  Are 
Same  as  Shown  in  Figure  B-l 


Dj  =  -  7  a/S 
n7  =  -  5  a/8 
n3  =  -  3  a/8 
n4  =  -  a/3 

^5  =  '  n4 

n6  =  -  n3 

n7  =  -  n2 


n8  =  "  nl 


(B-27 


The  present  calculation  is  restricted  to  the  case  when: 


APPENDIX  C 


GLIDE  SLOPE  SCATTERING  FROM  HORIZONTAL 
CYLINDRICAL  SURFACES 

The  computer  program  developed  by  the  Ohio  University  (O.U.) 
group**  to  treat  glide  slope  signal  scattering  by  aircraft  fuselages 
unfortunately  requires  excessively  long  running  times  for  its 
execution  even  on  the  most  rapid  digital  computers.  In  an  effort 
to  reduce  the  execution  time  of  this  program,  the  TSC  group  has 
been  investigating  approximate  techniques  for  evaluating  the  various 
integrals  which  appear  in  the  O.U. ,  formulation  of  the  fuselage 
scattering  problem.  Some  progress  has  been  made  in  this  area  and 
some  of  our  preliminary  results  will  now  be  presented. 

The  O.U.  researchers  simulate  the  fuselage  of  an  aircraft  by 
a  circular  cylinder.  In  formulating  the  problem  of  glide  slope 
signal  scattering  from  such  a  cylinder,  they  are  required  to 
evaluate  a  series  of  integrals  of  the  following  general  form: 


J  -  /  DM)  eitaf<«d*  , 


(c-i) 


where  k  =  2 a/A,  X  is  the  glide  slope  wavelength  (-  5  ft.j,  a  is 
the  radius  of  the  cylinder,  and  $  is  the  angular,  circular,  cylin¬ 
drical  coordinate.  The  numerical  evaluation  of  thest  :  tegrals 
is  very  time  consuming.  Our  group  has  found,  however,  that,  in 
some  cases,  they  can  be  evaluated  approximately  with  considerable 
accuracy  using  the  method  of  stationary  phase.  The  savings  in 
time  resulting  from  the  use  of  the  method  of  stationary  phase  as 
opposed  to  direct  numerical  integration  are  considerable. 

The  method  of  stationary  phase  is  a  technique  for  obtaining 
asymptotic  series  representations  for  integrals  of  the  type  given 
in  Equation  (C-1)  when  (ka)  is  a  large  parameter.  These  series 
expansions  arc  in  powers  of  (ka)’*.  The  technique  is  based  upon 
the  fact  that  when  (ka)  is  large,  the  principal  contributions  to 
the  integral  come  from  the  vicinities  of  the  end  points  of  the 
range  of  integration  (i.e.,  0  =  it/ 2  and  o  =  5n/2)  and  from  the 
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vicinities  of  those  points  $o  within  the  range  of  integration  at 
which  the  phase  function  f (<*>)  is  stationary: 


-  0 


(C-2) 


A  detailed  development  of  the  method  of  stationary  phase  can  be 

12 

found  in  the  monograph  Asymptotic  Expansions  by  A.  Erdelyi. 

Applying  Erdelyi 's  general  results  to  integrals  of  the  type  in 
Equation  (C-i),  we  obtained  the  following  approximate  representa¬ 
tion  which  applies  when  f($)  has  only  one  stationary  point  $>o  with¬ 
in  the  range  of  integration  and  when  f**($o)  <  0: 


J  *  A,  (ka) 


'1/2  +  A,  (ka)"1  +  A,  (ka)"3/2 


(C-3) 


where : 


si  ■  (tot)  1/2  Dt*o> 


i(kaf(^o)-a/4) 


a2  ■  1  [?%7tr  •“ . rwfr  e 


A.  - 
o 


D(tt/2)  rtikaf (»/2)  .  P(3«/2;  oikaf(3s/2)J 

S/2  (f<3)M2] 


A3,  pf.-s  ■)  iA 

4  UUoJ  2 


V  f(4)(V 


Jf"(*0)i5/2  6  lf"t 

1/2  f (3)c*j 


_ I 

6  M7/2J 


r=  /  \  1/2 

•  —  (rr-fyl) 


(C-4) 


(C-5) 


(C-6; 


If  f"(4>0)  is  positive,  the  factors  e*itt^  and 


which 

appear,  respectively,  in  Equation  (C-4)  and  (C-6)  should  be  replaced 
by  their  complex  conjugates.  Equations  (C-3),  (C-4).  (C-5),  and 
(C-6)  represent  the  leading  terras  in  the  asymptotic  series  expan¬ 
sion  of  the  integral  J  for  large  values  of  ka.  The  first  and  third 
terms  in  Equation  (C-3)  represent  contributions  to  J  from  the 
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& 
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vicinity  of  the  stationary  point  while  the  second  term  repre¬ 
sents  contributions  from  the  vicinities  of  the  end  points  -/ 2  and 
5a/2.  As  one  might  infer  from  the  complexity  of  the  coefficient 
A,  in  Equation  (C-6),  the  generation  of  any  additional  terms  in 
the  asymptotic  expansion  of  J  is  a  Herculean  task-  Nevertheless, 
we  have  been  able  to  obtain  reasonably  good  results  using  just  the 
three  terms  shown  in  Equation  (C-3).  As  an  example,  consider  the 
case  in  which  D(i»)  and  £(<P)  are  defined  as  follows: 


0($)  =  sin  -  HA  -  a  sin  ^sin 


kaz7 


sin  0^  , 


f($)  = 


C*Vha> 


sin  $ 


+  cos  3) 


cos  $  + 


aA‘ 

2p- 


cos“  $  - 


(C-7) 


(C-8) 


The  various  parameters  appearing  in  Equations  (C-7)  and  (C-8)  are 
defined  below: 


A 


=  height  of  the  axis  of  the  cylinder  above  the  ground 
=  hei-ht  of  the  transmitting  antenna  above  the  ground 


z2  =  height  of  the  receiving  antenna  above  the  ground 

8  =  aspect  angle  between  the  normal  to  the  surface  of  the 

cylinder  and  the  position  vector  of  the  receiver  rela¬ 
tive  to  the  cylinder 

r  =  distance  from  the  cylinder  to  the  receiver 

p  -  distance  from  the  transmitter  to  the  center  of  the 
cylinder 

A  =  p  multiplied  by  the  cosine  of  the  angle  of  incidence- 

let  Jj,  denote  the  value  of  the  integral  (C-i)  with  D  and  f 
defined  as  in  Equations  (C-7)  and  (C-8)  obtained  by  direct  numeri¬ 
cal  integration  using  Simpson's  rule  and  let  denote  the  cor¬ 
responding  asymptotic  approximation  to  the  integral  obtained  using 
the  series  expansion  (C-3)-  In  Table  C-l,  values  of  and  arc 
compared  for  several  values  of  the  aspect  angle  S  and  for  ka  =  21. 
The  running  times  required  on  the  PDP-iO  computer  to  obtain  each 
number  is  also  indicated. 
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TABLE  C-l.  NUMERICAL  INTEGRATION  AND  SERIES  EXPANSION  COMPARISONS 


8  =  80: 

JN 

=  - 

.783 

-  i(.824): 

time 

1320 

ms 

JA 

=  - 

.748 

-  i(.762) : 

time 

11.7 

ms 

3  =  60: 

JN 

=  - 

.531 

-  i(.801): 

time 

1360 

ms 

JA 

=  - 

.535 

-  i ( . 776) : 

time 

10 

ms 

3  =  40: 

JN 

=  - 

.337 

-  i(.480) : 

time 

1300 

ms 

JA 

=  - 

.329 

-  i(.468) : 

time 

10 

ms 

8  =  20: 

JN 

= 

.0868 

-  i(.S81) : 

time 

1400 

ms 

JA 

= 

.0867 

-  i(.594) : 

time 

11.7 

ms 

As  can  be  seen  fro®  the  table,  the  asymptotic  approximation 
yields  fairly  good  results  over  a  wide  range  of  aspect  angles.  The 
savings  in  time  which  result  from  the  use  of  this  approximate 
technique  are  clearly  enormous.  Much  more  work  will  be  required, 
however,  before  the  technique  is  fully  perfected. 
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APPENDIX  D 

REFLECTION  FROM  UNEVEN  GROUND  IN 
FRONT  OF  A  LOCALIZER  ANTENNA  ARRAY 


In  this  Appendix  we  develop  a  means  of  approximately  treating 
localizer  scattering  when  uneven  ground  exists  in  front  of  the 
localizer  such  as  a  hump  in  the  runway. 

When  a  localizer  antenna  is  situated  a  distance,  H,  above  a 
perfectly  conducting  flat  ground  plane  of  infinite  extent,  the 
electric  field  at  a  remote  observation  point  can  be  expressed  as 
(Reference  1,  Equation  (2.58)): 


=  u  Eo£($) 


k2F/Dp 


(D~l) 


where  u  is  the  unit  vector  giving  the  direction  of  polarization  of 
the  field;  f($)  is  the  horizontal  antenna  gain;  R  is  the  distance 
from  the  antenna  to  the  observation  point;  k  is  2«  over  the  wave¬ 
length;  Dp  is  the  projection  of  R  on  the  ground  plane;  H  is  the 
distance  of  the  antenna  from  the  ground  plane;  Eq  is  an  electric 
field  amplitude  and  z*  is  the  distance  of  the  observation  point 
from  the  ground  plane.  The  above  formula  for  is  restricted  to 
the  case  when  the  elevation  angle  of  the  ray  to  the  observation 
point  is  small. 


In  actual  airports,  however,  the  reflecting  ground  directly 
in  front  of  the  localizer  antenna  may  be  uneven.  Figure  D-l  shows, 
for  example,  the  lateral  profile  of  a  portion  of  the  ground  in 
front  of  the  localizer  antenna  at  the  Svracuse-Mancock  Airport. 


£ 
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In  order  to  treat  this  and  similar  cases  approximately  we 
consider  a  simple  model  where  the  uneven  foreground  consists  of  a 
number  of  connected  flat  surfaces  having  no  height  variation  in 
the  direction  perpendicular  to  the  vertical  plane  through  the 
runway  center  line.  The  end  of  the  runway  near  the  localizer  is 
taken  as  the  origin  of  the.-  coordinate  system  x-z.  Figure  D-2 
shows  a  diagram  of  the  simple  model  we  are  considering.  How  well 
this  simple  geometrical  model  works,  depends,  of  course,  on  the 
complexity  of  the  contours  of  the  actual  ground  being  considered. 
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Figure  D-l.  Lateral  Profile  Showing  a  Portion  of  the 

Ground  in  Front  of  the  Localizer  Antenna  at 
the  Syracuse-Hancock  Airport 
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Diagraa  Showing  a  Model  of  the  Reflecting  Ground 
in  Front  of  the  Localizer  Antenna.  The  End  of 
the  Runway  Near  the  Localizer  Is  the  Origin  of 
the  Coordinate  Syszes 
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In  order  to  still  be  able  to  use  Equation  CD-I)  for  this  uneven 
ground  plane  case,  we  consider  the  following  model.  We  assume 
geometric  optical  rays  emanating  from  the  localizer  antenno  and 
being  reflected  from  a  flat  ground  slab  segment.  The  rays  reflect¬ 
ing  from  the  two  ends  of  the  slab  define  a  reflection  region  on 
the  glide  path  as  shown  in  Figure  D-3.  Within  this  region  we  use, 
in  this  approximation.  Equation  D-l  with  H  now  the  perpendicular 
distance  of  the  antenna  from  the  extended  plane  of  the  ground 
slab;  and  z*  now  the  perpendicular  distance  of  the  observation 
point  on  the  glide  path  from  the  extended  plane  of  the  slab. 

Given  the  coordinates  of  the  end  points  of  the  slab  (x- .Sj) 
and  (xj, z-)  where  j  =  i+1  (see  Figure  D-3)  we  obtain  the  limits, 
XK-n  and  x-,^,  of  the  radiated  zone  on  the  glide  path: 
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xain  “  XGS  +  l2 
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Figure  D-4.  Diagram  for  Calculating  tfee  Quantities  H  and 
s’  Khcn  the  Observation  Point  Lies  on  the 
Overlapping  or  Gap  Regions 
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For  the  b  considered  above,  the  projection  of  the  localizer 
antenna  onto  the  extension  of  the  slab  has  a  perpendicular  dis¬ 
tance  given  by  (see  Figure  D-4): 


"  =  <xrxa>  (x]  -  :;) 


*  <V: 


cos  0- 


The  projection  of  an  observation  point  (x,z)  located  on  the  glide 
path  onto  the  extension  of  the  slab  has  a  perpendicular  distance 
given  by  (see  Figure  D-4): 


Az  cos  6. 


where: 


Az  =  (x-x-e)  Tan  -  (x 


^j'^up) 


Since  successive  slabs  make  different  angles  with  respect  to 
a  horizontal  ground  plane,  the  radiated  regions  of  successive  slabs 
may  overlap  or  there  may  be  gaps  between  them.  For  these  gaps  or 
overlapping  zones,  we  use  another  approximation  scheme  to  calculate 
the  quantities  H  and  z’.  Figure  D-4  shows  a  diagram  for  this 
method.  The  point  P  =  (Xp.Zp)  is  where  two  adjacent  slabs  meet. 

We  draw,  then,  a  line,  L2>  from  the  antenna  to  p  and  a  line,  Lj, 
from  p  to  (x,z)  a  point  on  the  glide  path  in  the  zone  of  the  gap 
or  overlap.  We  would  like  to  construct  an  imaginary  reflecting 
plane  for  these  rays.  Hence,  on  a  plane  tangent  to  the  rays  we 
draw  a  line,  pB,  through  the  point,  p,  bisecting  the  angle  between 
L2  and  Lj.  Then  the  plane  through  p  perpendicular  to  the  line  pB 
is  the  imaginary  reflecting  plane  required  for  the  rays.  For  the 
given  point  (x,z)  the  projections  H  and  z’  on  this  reflection  plane 
can  be  calculated  as  follows: 


x  -  xf 


hQ  =  lQ  Tan  <}>, 
z  “  ho  *  zup’ 
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This  set  of  equations  can  be  used  to  calculate  the  correspond¬ 
ing  H’s  and  z''s  for  a  sequence  of  points  in  the  gap  or  overlapping 
region.  We  note  that  for  a  convex  reflecting  bump  there  is  a  gap 
and  for  a  concave  reflecting  depression  there  is  an  overlap.  Also 
we  note  that  as  the  point  (x,z)  moves  to  the  boundary  points  of 
the  gap  region,  the  imaginary  plane  will  then  coincide  with  the 
actual  ground  planes  meeting  at  p. 

In  summary,  in  order  to  treat  approximately  the  scattering 
problem  when  non-flat  terrain  exists  in  front  of  the  localizer, 
we  begin  with  a  model  of  a  series  of  connected  slabs.  From  each 
slab  we  obtain  a  geometric  optics  reflection  region  on  the  glide 
path.  Then  we  find  regions  of  gaps  and  overlaps  for  which  we  use 
the  second  method  described.  For  the  remaining  regions  the  first 
method  described  is  used. 


1C  3 


APPENDIX  E 


THEORETICAL  LINEAR  ARRAY  PATTERNS 

The  CDI  pattern  that  is  measured  in  orbital  flight  around  a 
localizer  installation  can  be  quite  different  from  the  ideal  pat¬ 
tern  expected  on  the  basis  of  the  signal  distribution  to  the  array 
and  the  known  element  patterns.  In  many  instances,  this  may  be 
only  partially  explained  by  the  influences  of  the  scattering 
environment.  An  important  additional  factor  can  be  mutual  couplings 
between  the  array  elements,  the  magnitudes  of  which  can  be  specific 
to  the  particular  localizer  installation.  The  consequence  of  this 
is  that  spurious  signal  currents  are  induced  in  individual  elements 
which  alter  the  total  currents  from  the  values  implied  by  power 
dividing  networks,  phase  shifters,  etc.  Instrumentation  is 
available  to  measure  actual  element  current  magnitudes  and  phases. 
Given  such  measurements,  it  is  possible  to  calculate  the  signal 
patterns  that  would  be  radiated  by  the  array  situated  on  a  flat 
ground  plane.  Routines  for  accomplishing  these  calculations  for 
an  arbitrary  linear  array  localizer  have  been  incorporated  in  the 
ILSLOC  program^  and  are  called  into  use  by  designating  the  variable 
"MODE"  as  7  or  8.  Detailed  input  data  specifications  are  given  in 
the  ILSLOC  user's  manual.3 
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